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EXECUTIVE  SUMMARY/INTRODUCTION 

As  a  first  step,  PCKA  and  NSWC  met  and  mutually  agreed  upon  an  example  system  to  be 
simulated  for  the  purpose  of  investigating  high-power  machinery  systems.  This  system,  which  is 
shown  later  in  Fig.  1,  is  similar  to  the  proposed  integrated  power  system  (IPS)  in  that,  the  salient 
features  of  the  IPS  are  included.  Briefly,  the  sample  IPS  considered  consists  of  a  3-phase 
synchronous  machine,  ac  distribution,  a  ship  service  power  supply,  a  15-phase,  19  MW 
propulsion  drive  and  a  pulsed  power  load.  This  system  was  selected  to  illustrate  the  utility  and 
facility  of  simulation  in  the  analysis  and  design  of  high-power  machinery  systems. 


Ship  Service  (DC) 
1  MW 


Figure  1 .  System  studied. 


In  the  original  SBIR  Phase  I  proposal,  PCKA  suggested  that  the  selected  system  be 
simulated  in  both  detailed  and  reduced-order  models  and  that  the  concept  of  waveform 
reconstruction  be  used  whenever  appropriate.  This  technique  (waveform  reconstruction),  which 
was  introduced  by  PCKA  in  a  previous  SBIR  Phase  II,  provides  a  means  to  accelerate  the  speed 
of  computation  by  extracting  information  from  the  reduced-order  model  to  reconstruct  the 
instantaneous  waveforms  of  the  system  variables.  However,  during  the  interim  between 
submitting  the  proposal  and  the  contract  award,  PCKA  personnel  were  successful  in  developing  a 
simulation  technique  which  markedly  increased  the  speed  of  computation  of  the  detailed  system 
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simulation  [1].  As  this  technique  developed,  it  became  clear  that  waveform  reconstmction  of  the 
system  variables  from  the  reduced-order  simulation  was  less  important.  Therefore,  this  approach 
was  not  pursued  even  though  this  technique  could  be  applied  using  the  reduced-order  models 
developed  in  this  SBIR.  Instead,  attention  was  focused  on  demonstrating  the  facility  of  a 
detailed  and  reduced-order  system  simulation  for  analysis  and  design  purposes  of  high-power 
systems. 

In  this  SBIR  Phase  I,  the  detailed  and  reduced-order  models  of  the  example  IPS  were 
developed.  Example  computer  studies  are  given  to  illustrate  the  system  performance  during 
various  large-disturbances  including  the  operation  of  the  pulsed  power  weapon.  For  purposes  of 
validating  the  reduced-order  modeling  techniques,  duplicate  studies  were  performed  using  the 
detailed  and  reduced-order  system  simulations.  These  studies  are  presented  in  this  report  along 
with  an  analytical  description  of  the  simulation  of  each  system  component.  In  some  cases, 
however,  references  are  cited  rather  than  repeating  previous  derivations.  The  digital  code  is  also 
given  to  the  Program  Manager,  Mr.  Henry  Robey  of  NSWC  Annapolis,  so  that  these  simulations 
(detailed  and  reduced-order)  can  be  used  by  the  Navy  and/or  their  contractors. 

There  are  additional  aspects  of  this  work  which  are  considered  innovative  with  applications 
well  beyond  the  proposed  IPS.  In  particular,  the  analysis  of  the  15-phase  induction  machine  is 
new  and  the  simulation  developed  should  have  application  to  other  multiphase  induction  or 
synchronous  machines  [Appendix  B].  Another  aspect  of  this  work  deals  with  system  stability.  It 
was  found  that  the  field-oriented  control  (vector  control)  associated  with  the  proposed  15-phase 
propulsion  induction  motor  drive,  exhibited  negative-impedance  instability  [Appendix  C].  This 
was  verified  by  implementing  in  hardware  a  lower  power,  3-phase  induction  motor  with  the  same 
type  of  vector  control  [Appendix  A].  It  was  shown  that  this  system  also  inhibited  negative- 
impedance  instability.  In  other  words,  this  verified  that,  when  placed  in  a  system,  vector 
controlled  machines  could  promote  negative-impedance  system  instability  in  a  manner  similar  to 
that  of  constant  power  loads.  To  our  knowledge,  this  has  not  been  demonstrated  previously. 
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Perhaps  more  important  than  demonstrating  that  negative-impedance  system  instability 
could  occur  due  to  vector  controlled  machines,  is  the  fact  that  PCKA  personnel  (S.  D.  Sudhoff) 
developed  a  control  technique  which  prevents  this  instability  [Appendix  A].  In  particular,  by 
incorporating  an  additional  control  within  the  vector  control,  system  stability  can  be  maintained 
with  insignificant  degradation  of  the  performance  of  the  vector  controlled  machine.  It  is  of 
utmost  importance  to  emphasis  that  the  negative-impedance  feature  of  the  vector  control  was 
discovered  through  the  detailed  simulation  of  the  proposed  IPS.  Moreover,  the  technique  to 
eliminate  this  instability  was  also  developed  through  simulation  before  implementing  the 
hardware  prototype.  These  facts  alone  attest  to  the  success  of  this  SBIR  Phase  I  research  by 
demonstrating  the  value  of  computationally  efficient  detailed  models  and  nonlinear  reduced- 
order  models  derived  from  and  validated  against  these  detailed  models.  This  work  sets  the  stage 
for  Phase  II  wherein  an  ACSL-based  simulation  library  for  shipboard  electric  power/drive 
systems,  specifically  the  IPS,  will  be  developed  and  demonstrated.  In  Phase  III  this  library 
would  be  made  available  by  the  Navy  in  the  advanced  surface  and  machinery  program  (ASMP) 
and  to  Navy  contractors  for  the  IPS  full  scale  engineered  development  (FSED)  system 
integration  efforts. 


SYSTEM  DESCRIPTION  AND  MODEL  DEVELOPMENT 

The  system  considered  in  this  smdy  is  illustrated  in  Fig.  1 .  The  turbine  is  the  prime  mover 
for  a  3-phase,  60  Hz,  21  MW,  4160  V  (line-to-line)  synchronous  machine.  The  voltage  is 
controlled  by  an  IEEE  Type  2  excitation  system.  The  synchronous  machine  is  connected  to  a  3- 
phase  4160  V  ac  distribution  bus.  The  principal  load  on  the  bus  is  the  19  MW  ship  propulsion 
motor  drive.  In  addition  to  the  propulsion  load,  there  is  a  1  MW  ship  service  power  supply 
consisting  of  a  12-pulse  controlled  rectifier,  as  well  as  a  pulsed  power  load.  A  harmonic  filter  is 
present  which  reduces  the  harmonics  caused  by  the  electronic  switching  in  the  motor  controllers 
and  ship  service  power  supply.  The  bus  capacitance  shown  represents  the  stray  bus  capacitances. 
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A  detailed  description  of  each  of  these  components  follows,  along  with  a  brief  discussion  of  the 
average-value  modeling  strategy  as  appropriate. 

Prime  Mover 

Since  the  focus  of  this  effort  was  the  representation  of  the  electrical  dynamics,  the  prime 
mover  was  assumed  to  be  ideal  and  therefore  the  synchronous  machine  speed  was  assumed  to  be 
constant  and  the  frequency  60  Hz. 

Synchronous  Machine 

The  synchronous  machine  is  a  2-pole,  21  MW  (26.25  MV A),  4160  V  (line-to-line),  salient 
pole  machine.  The  parameters  are  given  in  Table  1.  The  modeling  of  a  synchronous  machine 
has  been  reported  previously  by  PCKA  and  therefore  it  is  not  repeated  here  [2]. 


Table  1.  Synchronous  machine  parameters 

P  =  2 

Tj  =  1.27mQ 

4=391  /zH 

4^  =  2.51/zH 

r;=5.26mn 

4^  =  157  /zH 

4,=2.79azH 

r^=4.74ma 

4,  =  69.8;zH 

r'^j  =  40\mQ. 

4,  =227  /ZH 

M„=  3.73.6  H-' 

=360.5  H-' 

Tj.  =  1.90  H-' 

Excitation  System 

The  synchronous  machine  excitation  system  /  voltage  regulator  is  an  IEEE  Type  2.  The 
block  diagram  is  given  in  Fig.  2  and  the  parameters  are  listed  in  Table  2.  These  are  reasonable 
parameters  based  on  the  size  of  the  machine.  Unfortunately,  the  parameters  provided  by  the 
manufacturer,  appear  to  be  incorrect. 
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Figure  2.  IEEE  Type  2  exciter  representation. 


Table  2.  Exciter  parameter 

=  4160  V 

^RMAX  “  7-3  P-U. 

S^„=0.50 

=  0.8  s 

^E]oo  ~  0‘S6 

=  400 

Kp  =  0.03 

II 

b 

=  20  ms 

=  1  s 

=  0.8  s 

Tp,  =  1  s 

Induction  Motor 

The  induction  motor  which  is  used  for  ship  propulsion  is  a  150  rpm,  19  MW  unit.  The 
steady-state  equivalent  circuit  is  shown  in  Fig.  3  and  the  parameters  for  the  five  3-phase  set 
model  are  given  in  Table  3.  Modeling  of  the  15-phase  induction  motor  was  a  major  technical 
challenge.  Since  this  work  may  be  of  interest  to  a  relatively  large  audience,  a  description  of  the 
modeling  is  documented  in  Appendix  B. 
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Figure  3.  Per-phase  equivalent  circuit  of  15-phase  induction  machine. 


Table  3.  Induction  motor  parameters. 

4=6.79  mH 

4  =  4.57  mH 

r,  =  122  mQ 

4  =58.1mH 

«,r=  0.373 

P  =  12 

Motor  Controller 

The  induction  motor  power  converter,  which  converts  4160  V,  3-phase  ac  to  variable 
amplitude,  variable  frequency,  15-phase  ac  is  shown  in  Fig.  4.  The  converter  consists  of  three 
uncontrolled  rectifiers  connected  to  three  dc  links.  Each  dc  link  is  connected  to  five  H-bridge 
inverters  (Fig.  5)  and  each  H-bridge  supplies  one  phase  of  the  induction  motor  (Fig.  4).  The 
parameters  are  given  in  Table  4. 


Table  4.  Induction  motor  power  converter  parameters. 

Symbol 

Description 

Value 

Effective  thyristor  on  state  resistance 

10  m£2 

Effective  thyristor  off  state  resistance 

100  Q 

Load  commutated  converter  commutating  inductance 

10  fiR 

^dc 

DC  link  inductance 

3mH 

DC  link  resistance  (associated  with  4,-) 

10  iiQ. 

DC  link  capacitance 

10  nF 

Effective  diode  voltage  drop 

5  V 

Effective  thyristor  voltage  drop 

8  V 
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Figure  4.  Induction  motor  power  converter. 


Figure  5.  H-bridge. 
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The  most  difficult  technical  challenge  of  this  project  was  the  development  of  the 
appropriate  algorithms  to  control  each  of  the  semiconductors  in  the  15  H-bridges.  According  to 
the  manufacturer,  the  overall  control  strategy  was  to  be  based  upon  field  oriented  control. 
Furthermore,  the  control  strategy  was  to  be  implemented  in  such  a  way  that  the  switching 
frequency  of  the  H-bridge  semiconductors  was  to  be  constant.  There  were  two  principal 
difficulties  in  implementing  such  a  control.  First,  since  the  machine  is  a  15-phase  design,  the 
interaction  between  the  motor  phases  complicated  the  design  of  the  current  control  loop. 
Secondly,  and  most  importantly,  it  was  discovered  that  if  a  standard  field  oriented  control  was 
used  the  induction  motor  propulsion  drive  would  cause  the  system  to  become  unstable. 
Therefore,  a  new  concept  which  is  denoted  "link  stabilizing  field  oriented  control"  was 
developed.  This  control  acts  much  the  same  as  a  standard  field  oriented  control  with  the 
exception  that  it  does  not  cause  the  drive  to  exhibit  the  negative  impedance  instability  as  does  the 
standard  field-oriented  control.  This  new  concept  in  induction  motor  control  has  ramification’s 
far  beyond  this  system.  It  is  described  in  detail  in  Appendices  B  and  C. 

Ship  Service  Power  Supply 

The  ship  service  power  supply  is  depicted  in  Fig.  6.  Therein,  i^p,  and  denote  the 

x'th  phase  currents  of  the  primary,  delta-connected  secondary,  and  wye-connected  secondary 
transformers,  respectively;  the  variables  v^^p,  and  designate  the  a  to  b  primary, 

delta-connected  secondary,  and  wye-connected  secondary  transformer  voltages,  respectively. 
The  rectifier  voltage  and  current  are  denoted  v^  and  i^.  The  transformer  output  voltage  and 

current  are  denoted  vdc  and  idc-  The  individual  SCR’s  of  the  upper  load  commutated  converter 
are  denoted  D1-D6  and  the  SCR’s  of  the  lower  converter  are  W1-W6.  The  rectifiers  are 
connected  in  series  and  an  LC  filter  is  connected  between  the  rectifier  output  and  the  load. 
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V  ^dc 


Figure  6.  Power  supply. 


The  supervisory  control  algorithm  for  the  power  supply  is  depicted  in  Fig.  7.  Therein,  the 
PI  voltage  control  is  the  principal  control  loop  with  a  current  control  to  prevent  over  corrects. 
The  thyristor  gate  signals  are  fired  based  on  the  positive  zero  crossing  (PZCD  in  Fig.  7)  of  the  a 
to  b  delta  connected  secondary  voltage.  In  particular,  the  individual  thyristor  valves  are  fired  in 
accordance  with  Table  5. 


Table  5.  Thyristor  turn  on  times 

SCR 

Tum-on  Time 

SCR 

Tum-on  Time 

D1 

a/  (0^ 

W1 

{a  +  \\K !  €) !  0)^ 

D2 

(a  +  f)/tD, 

W2 

{a  +  Kl€)l  (o^ 

D3 

(«  +  ¥)/ 6). 

W3 

{a  +  3Tcl  6)1  (0^ 

D4 

(a  +  ¥)/®. 

W4 

{a  +  5Kl  6)1  co^ 

D5 

[a  +  ^)l(o^ 

W5 

{a  +  lKl  6)1  (0^ 

D6 

{a  +  ^)l  0)^ 

W6 

{a-r9nl  6)1  (0^ 
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Figure  7.  Power  supply  control  algorithm. 

Parameters  of  the  ship  service  power  supply  are  listed  in  Table  6;  therein  the  transformer 
parameters  are  those  of  the  simplified  equivalent  circuit  shown  in  Fig.  8.  Note  that  the  primary 
leakage  inductance  is  negative;  this  results  from  the  fact  that  the  primary  leakage  in  Fig.  8  is 
actually  an  effective  value  based  on  several  terms,  not  because  any  physical  leakage  is  claimed  to 
be  negative  [2]. 
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Table  6.  Power  supply  parameters 


Description 


Lu  Secondary  leakage  inductance 


Value 


28.2  mH 


325  mQ 


Primary  effective  turns  /  Delta-connected  secondary 
effective  turns  (ratio  to  wye -connected  secondary  is 
selected  to  give  the  same  line-to-line  voltage) 

9.8 

4 

Primary  leakage  inductance 

-9.69  mH 

Primary  resistance 

274  mQ 

^dc 

Filter  inductance 

600  jUH 

^dc 

Filter  inductance  resistance 

1  mQ 

Cdc 

Filter  capacitance 

3200  AtF 

Thyristor  off  state  resistance 

100  a 

Thyristor  on  state  resistance 

1  a 

K,- 

Voltage  controller  proportional  gain 

0.01 

Kip 

Voltage  controller  integral  gain 

20/s 

limv 

Voltage  regulator  limit 

lOkV 

K, 

Current  regulator  proportional  gain 

0.01  V/A 

Kd 

Current  regulator  integral  gain 

0  V/As 

limi 

Current  regulator  limit 

lOkA 

amin 

Minimum  phase  delay 

0 

^max 

Maximum  phase  delay 

3.14  radians 

V 

max 

Maximum  power  supply  output  voltage 

1147  V 

T 

Line-to-line  voltage  filter  time  constant 

4.594  ms 
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Figure  8.  Power  supply  transformer  simplified  equivalent  circuit. 


The  detailed  computer  model  of  the  power  supply  was  developed  under  previous  Navy 
contracts  [2].  The  average- value  model  used  in  the  studies,  reported  herein,  was  developed 
under  this  SBIR.  Therefore,  it  is  appropriate  to  set  forth  briefly  the  modeling  strategy.  The 
average-value  model  of  the  power  supply  is  based  upon  the  average-value  model  of  a  voltage 
behind  inductance  /  controlled  rectifier  model.  Although  the  system  is  a  12-pulse  rectifier,  it  is 
convenient  to  consider  the  operation  of  a  6-puIse  load  commutated  converter,  as  depicted  in  Fig. 
9,  before  proceeding  to  the  12-pulse  rectifier.  Since  an  average-value  model  is  desired,  the 
known  quantities  will  be  considered  to  be  the  source  voltage  (in  qd  form)  and  the  dc  output 
voltage  V  The  output  of  the  model  will  be  the  fast-average  q-  and  d-axis  source  currents 

and  the  fast-average  of  the  rectifier  current  [3].  For  this  purpose,  it  is  convenient  to  assume  that 
the  source  voltage  may  be  expressed  in  terms  of  q-  and  d-  axis  voltages  as 


(1 
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where  E  is  the  rms  magnitude  of  the  source  voltage.  The  d-axis  voltage  is  zero;  this  can  be 

obtained  by  appropriate  selection  of  the  reference  frame.  The  reference  frame  is  the  source 
reference  frame  indicated  by  a  superscript  g  and  an  angular  displacement  9^ . 

In  order  to  obtain  the  fast-average  rectifier  voltage  and  the  q-  and  d-axis  currents,  it  is 
convenient  to  note  that  these  variables  are  periodic  in  7t/3  radians  of  0^.  Thus,  the  average  may 
be  taken  over  any  Tt/3  interval  of  9^.  It  is  convenient  to  consider  the  an  interval  which  begins 

when  a  certain  valve  turns  on  and  continues  until  the  next  valve  is  turned  on.  If  the  interval  when 
valve  3  begins  to  conduct  until  valve  4  begins  to  conduct  is  selected 

2;r 

— +ot 

=  -  j  (3) 

— +a 
i 

In  (3),  a  is  the  firing  angle  relative  to  the  source  reference  frame.  In  particular,  a  is  defined  such 
that  valve  3  fires  when 

0=-  +  a  (4) 

^  3 

The  average  indicated  in  (3)  may  be  carried  out  by  noting  from  Fig.  9 


(5) 
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V,  =  V,  +  L 
fa  %  .c 


V  =  V  +  L  — ^ 

«  c,  ac 


(6) 

(7) 


Substitution  of  (5)-(7)  into  (3)  yields 

=|izr  +^4c^y,(4 -4)li 

^  i 


(8) 


Equation  (8)  may  be  simplified  by  transforming  (1)  and  (2)  back  to  abc- variables,  whereupon 


=  V2  Ecos^ 


'^hg 


=  V2Ecos(0^-^) 


v^,^,  =  V2  Ecos(0^  +^) 


(9) 

(10) 

(11) 


Substitution  of  (9)  -  (1 1)  into  (8)  yields 

_  3^|3 


llZ 


V.  = 


V2  E  cos  a  +  — 4^6)  (4-4)1;,^ 
K  K  ^  - 


— +a 
3 


(12) 


Further  simplification  may  be  obtained  by  observing  that,  prior  to  the  instant  valve  3  begins  to 
conduct,  only  valves  1  and  2  are  conducting;  therefore 


K  0  If 


(13) 


In  (13),  dc  current  ripple  is  neglected  whereupon  the  instantaneous  rectifier  current  is  equal  to  it's 
fast-average  value.  Similarly,  at  the  instant  prior  to  valve  4  conducting  only  valves  2  and  3  are 
on,  therefore 

|4faZ  =  [0  - 4 - I  +  (^4) 


15 


In  (14),  A/,  represents  the  change  in  the  fast  average  of  the  rectifier  current  over  the  interval.  It 
follows  from  this  definition  that  the  derivative  of  the  fast-average  of  rectifier  current  may  be 
approximated  as 


di^ 

It 


(15) 


Substitution  of  (13)-(15)  into  (12)  yields  the  desired  expression  for  the  fast-average  rectifier 
voltage, 


_  3  -  di 

^  V2  £cos  a  -  - -  24, 


(16) 


It  is  assumed,  that  the  rectifier  is  connected  to  a  dc-link  inductor  with  inductance  and 
resistance  rjc-  This  output  voltage  is  denoted  (Fig-  9).  The  voltage  equation  becomes 

di. 


Vr  =  >'Jr  +  Lj,-^Vj, 


(17) 


The  fast-average  of  (17)  is 


_  -  di 

Vr=>-jJr+Luc-^^dc 


(18) 


Combining  (16)  with  (18)  yields 

,j  —  V2  £  cos  a  -  +  7  4cCyJv  - 

^-_rc _ _ _  (19) 

dt  4^+24^ 

It  is  necessary  to  calculate  the  average  q-  and  d-axis  source  currents.  In  order  to  formulate 
the  expressions  for  the  average  q-  and  d-axis  currents,  it  is  assumed  that  the  rectifier  current  is 
constant  throughout  the  interval  and  equal  to  its  fast-average  value.  It  is  convenient  to  divide  the 
interval  into  two  parts  in  order  to  compute  the  average  current;  the  commutation  interval  dunng 
which  the  current  is  transferred  from  valve  1  to  valve  3,  and  the  conduction  interval  during  which 
only  valves  2  and  3  are  conducting.  During  the  commutation  interval,  valves  1,  2,  and  3  are 
conducting.  Therefore,  the  current  into  the  ac  source  must  be  of  the  form 
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Lhcs  4  ^as  ^r] 


V  —  Vj,  =  0 
as  bs  ^ 


Substitution  of  (5),  (6),  and  (20)  into  (21)  and  solving  for  the  a-phase  current  yields 

cos(a)-cos(0^--|)  (22) 

The  commutation  sub-interval  ends  when  the  current  in  valve  1,  which  is  the  a-phase 
current,  becomes  zero.  The  angle  from  the  time  valve  3  is  turned  on  and  valve  1  is  turned  off  is 
known  as  the  commutation  angle  u.  It  can  be  found  by  setting  (22)  equal  to  zero.  In  particular. 


u  = -a  +  arccos  cos  a  - 


L  CO  I 


After  commutation,  the  a-phase  current  remains  at  zero,  therefore 

U.=[0-i;tr  (24) 

Equations  (20),  (22),  and  (24)  specify  the  machine  currents  throughout  the  kII  interval  which 
begins  when  valve  3  is  fired. 

The  next  step  is  to  use  the  machine  currents  to  calculate  the  average  q-  and  d-  axis  currents, 
which  are  defined  as 


(2« 

Since  the  form  of  machine  currents  is  different  in  the  conduction  period  than  in  the  commutation 
period,  it  is  convenient  to  break  up  (25)  and  (26)  into  components  corresponding  to  those  two 


intervals.  In  particular. 


^qs,com  ^qs,cond 


^ds,com  ^dsxond 


where 


17 


(29) 

(30) 

(31) 

(32) 

The  commutation  component  of  the  current  may  be  found  by  substituting  (22)  into  (20), 
transforming  the  result  to  the  source  reference  frame  Q)g,  and  then  integrating  in  accordance  with 


(29)  and  (31) 


qsxom 


2V3  T 
= - K 

K 

3  a/3£ 


5it\  .  (  Stv 
sinl  u  +  a-  —  j  -  smi  a  -  — 


cos(a)[cos(M  +  a)  -  cos(a)]  + 


1  3  42E 
4  K 


[cos(2m)  —  cos(2a  +  2i/)] 


(33) 


2V3  r 


(Is,  com 


K 

3  V2£ 

K  4,0)^ 


5;r^  (  ‘S7c\ 

-cosi  M  +  a-  — l  +  cosi  a-  — 


cos(a)[sin(M  +  a)  -  sin(Q:)]  + 


1  3  V2£  r  . 


A  K  L  (0 

^ac  ^ 


rsin(2M)  -  sin(2a  +  2m)] - - 

K  L 


3  V2£  1 


,  (O  2 

■ac^g 


(34) 


Similarly,  the  conduction  component  of  the  average  currents  may  be  expressed  by  transforming 
(24)  to  the  source  reference  frame  followed  by  the  integration  as  set  forth  in  (30)  and  (32) 
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2V3  rf  .  f  iTt^  .  f  ,  ^57: 
- 1  sin  a  H - -  sin  a  +  m  h - 


i 


dsxond 


2^r[  1k\  I 

- 1  -cos  an - +COS 

TT  1  6  J  ' 


5;r^ 

Of  +  W  H - 

6  )\ 


(36) 


The  block  diagram  of  the  basic  rectifier  model  is  illustrated  in  Fig.  10.  Therein,  the  inputs 
are  the  firing  angle  a,  the  q-  and  d-axis  source  voltage  in  an  arbitrary  reference  frame  and  v“, 

and  the  dc  output  voltage  "^^^e  outputs  of  the  model  are  the  fast-average  of  the  rectifier 
current  and  the  fast-average  of  the  q-  and  d-  axis  currents  in  the  arbitrary  reference  frame 
and  Ij  .  The  calculation  in  Block  1  determines  the  peak  magnitude  of  the  source  voltage  V2F. 


Figure  10.  Diagram  of  rectifier  model  for  6-pulse  power  supply. 


The  calculation  in  Block  2  determines  the  angle  of  the  source  reference  frame  relative  to  the 
reference  frame  in  which  the  qd  voltages  are  inputs.  This  calculation  is  necessary  since  the  q- 
and  d-  axis  currents  are  calculated  in  the  source  reference  frame  (co^)  which  may  not  be  the 

reference  frame  in  which  the  voltages  are  input  to  the  rectifier  model.  With  the  magnitude  of  the 
source  voltage  and  the  dc-link  voltage,  the  time  derivative  of  the  rectifier  current  can  be 
calculated  (Block  3),  which  is  integrated  in  order  to  determine  the  fast-average  rectifier  current 
ij.  (Block  4).  With  the  rectifier  current,  the  magnitude  of  the  source  voltage,  and  the  firing  angle, 
the  commutation  angle  can  be  calculated  using  (23),  (Block  5).  Based  upon  the  rectifier  current. 
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the  commutation  angle,  and  the  firing  angle,  the  q-  and  d-axis  currents  are  calculated  in  the 
source  reference  frame  using  (27),  (28),  and  (33)  -  (36)  in  Block  6.  Finally  the  q-  and  d-  axis 
currents,  in  the  source  reference  frame,  are  transformed  to  the  arbitrary  reference  frame  in  Block 

7. 

The  rectifier  model  in  Fig.  10  is  for  a  6-pulse  rectifier  as  depicted  in  Fig.  9.  In  order  to 
modify  the  model  to  accommodate  the  12-pulse  rectifier,  is  should  be  noted  that  as  long  as  the 
commutation  angle  is  less  than  71/6,  the  two  6-pulse  rectifiers  will  not  interact.  In  other  words, 
each  rectifier  can  be  analyzed  separately,  ignoring  the  current  in  the  other  rectifier.  A  detailed 
explanation  of  the  interaction  between  6-pulse  units  can  be  found  in  [4].  The  changes  in  the 
structure  of  the  model  are,  (1)  the  output  voltage  is  doubled  since  there  are  two  rectifiers 

connected  in  series,  and  (2)  the  current  into  the  source  (the  primary  current)  is  also  doubled.  The 
model  is  illustrated  in  Fig.  11.  An  additional  modification  has  been  incorporated  in  order  to 
account  for  the  turns  ratio  of  the  transformer;  the  model  in  Fig.  1 1  is  referred  to  the  secondary  of 
the  transformer. 


Figure  11.  Diagram  of  rectifier  model  for  12-pulse  power  supply. 


In  order  to  complete  the  average-value  model,  it  is  necessary  to  integrate  the  rectifier  model 
depicted  in  Fig.  1 1  with  the  model  of  the  output  capacitor  shown  in  Fig.  6  and  to  add  the  model 
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of  the  controls  shown  in  Fig.  7.  It  is  clear,  that  the  capacitor  voltage  (in  average-value  form)  can 
be  calculated  from 


PVjc  = 


(37) 


An  interesting  aspect  of  the  modeling  is  representing  the  single-phase  filter,  the  positive 
zero  crossing  detector,  and  the  gate  sequencer  in  average-value  qd  form.  This  is  somewhat 
involved  since  the  firing  angle  is  determined  by  the  controller  (Fig.  7)  relative  to  the  filtered  a 


to  b  delta  connected  secondary  voltage  whereas  the  firing  angle  required  in  the  rectifier  model  is 
relative  to  the  primary  voltage.  The  firing  of  the  power  supply  thyristors  is  derived  from  the 
filtered  zero  crossing  of  the  aiob  delta  connected  secondary  voltage.  This  first-order  low-pass 
filter  has  a  time  constant  such  that  it  introduces  a  60°  phase  lag.  In  terms  of  abc  variables. 


For  the  purpose  of  representing  the  filter  as  an  average-value  model,  it  is  convenient  to  treat 
the  system  as  if  each  line-to-ground  voltage  is  filtered  individually  and  then  subtracting  the  b- 
phase  voltage  from  the  a-phase  voltage  in  order  to  determine  The  3-phase  set  of  filters 

may  be  represented  in  fast-average  qd  form  as 


dt  T 


where  and  vjj  are  the  q-  and  d-axis  delta  connected  secondary  voltages,  and  are 
the  corresponding  filtered  secondary  voltages,  and  (Og  is  the  speed  of  the  synchronous  reference 


frame. 

The  voltages  and  vjj  are  not  directly  available  from  the  model;  instead  they  are 
calculated  with  a  reduced-order  representation  of  the  circuit  depicted  in  Fig.  8.  In  particular,  Vqd 
and  vdc  may  be  estimated  as 
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^qd 


-l:+u 


•Ip  fdp 


/  n 


(41) 


Vdd  = 


—  r  i 
2 


(0. 


tip 


/  n 


(42) 


The  electrical  dynamics  are  neglected  in  (41)  and  (42);  however,  this  does  not  introduce  any 
noticeable  inaccuracy  since  the  time  constant  of  the  filter  is  relatively  large. 

The  filtered  a-  to  b-  phase  line-to-line  voltage  may  be  expressed 


v.w=V2v</Cos(^.  +  ^vj)  (43) 

where  Vqj  is  the  rms  amplitude  and 

<l>.d  =  ^^sH\d  -  J^dd)  (44) 

In  (44),  j  =  .  The  positive  zero  crossing  of  (43)  occurs  at 

e,  =  -k,-^  (45) 

The  gate  sequencer  will  turn  on  valve  3  of  the  delta-connected  rectifier  at  an  angle  of  a^+lK/3 
after  the  zero  crossing,  which  occurs  when  the  position  of  the  synchronous  reference  frame  is  at 

~  ^vd  (46) 

6 

Recall,  from  the  derivation  of  the  reduced-order  model,  valve  3  fires  at  an  angle  of  a^  +  K  !  3 
past  the  peak  value  of  the  a-phase  voltage,  which  is 

e,  =  |  +  a-^„,  (47) 

where 

=  angle(v^p -yvjp)  (48) 

Comparing  (47)  and  (46)  yields  the  expression  for  the  effective  firing  angle  relative  to  the  source 
voltage 

(^  =  cCc+Kp-kd-^  (49) 

A  block  diagram  of  the  average-value  model  of  the  gating  controls  is  shown  in  Fig.  12. 
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Figure  12.  Average-value  model  of  gate  control. 


Pulsed  Power  Weapon 

The  pulsed  power  weapon  considered  is  shown  in  Fig.  13.  It  consists  of  a  6-pulse  rectifier 
with  a  LC  pulse  forming  network  interfaced  to  the  system  by  a  wye-to-wye  transformer.  The 
rectifier  is  controlled  so  as  to  gradually  charge  the  capacitor  which  is  then  suddenly  discharged. 
The  control  of  the  pulsed  power  unit  is  divided  into  three  parts;  a  charge  profile  generator,  a 
firing  angle  control,  and  a  gate  sequencer. 


Figure  13.  Pulsed  power  charger. 


The  charging  profile  control  (Fig.  14)  generates  the  commanded  capacitor  voltage  as  a 
function  of  time.  As  shown  in  Fig.  14,  the  capacitor  voltage  increases  linearly  to  a  prescribed 
value  of  over  a  time  interval  r^/i-  The  voltage  command  then  remains  constant  for  This  brief 

period  is  used  to  allow  the  rectifier  to  finish  charging  the  capacitor  to  the  final  value.  After  this 
period  is  concluded,  the  firing  of  the  rectifier  is  disabled  and  gating  signals  cease.  The  period 

is  provided  to  allow  any  remaining  thyristor  valves  to  commutate  off  The  cycle  is  completed 
with  the  time  rj  during  which  the  capacitor  is  discharged. 
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* 


Figure  14.  Charging  profile  generator. 


The  firing  angle  control  varies  the  converter  firing  angle  so  that  the  capacitor  voltage  tracks 

the  desired  voltage  established  by  the  charging  profile  generator.  This  control  is  illustrated  in 
Fig.  15.  Therein,  the  commanded  capacitor  voltage  is  compared  to  the  actual  capacitor 

voltage  vcpp  and  the  difference  is  multiplied  by  a  gain  and  limited  to  give  the  current  command 
r  The  difference  between  the  commanded  rectifier  current  and  the  actual  current  is  then 

rpp' 

multiplied  by  a  proportional  gain  and  the  result  added  to  the  capacitor  voltage  to  obtain  the 
rectifier  voltage  command  The  capacitor  voltage  is  included  in  the  rectified  voltage 

command  since  it  decouples  the  capacitor  and  inductor.  Finally,  the  voltage  command  is  divided 
by  the  maximum  rectifier  voltage  Vrmax^  limited,  and  translated  into  a  commanded  firing  angle 

‘^cpp- 

'^cpp 

V 

cpp 

’VP 


Figure  15.  Firing  angle  control. 

The  gate  sequencer  is  illustrated  in  Fig.  16  and  is  considerably  different  from  that  used  for 
the  power  supply.  In  particular,  all  three  line-to-line  voltages  are  filtered,  and  the  firing  of 
individual  valves  is  timed  (via  the  delay  blocks)  from  either  the  positive  or  negative  zero 
crossings  of  line-to-line  voltages  as  indicated  in  Fig.  16.  Each  DELAY  block  represents  a  delay 
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of  radians  from  the  zero  crossing.  In  this  control  the  firing  angle  is  updated  6  times  every  cycle 
rather  than  once  a  cycle  as  in  the  case  of  the  power  supply.  Parameters  of  the  charging  profile 
generator,  the  firing  angle  control,  and  the  gate  sequence  are  given  in  Table  7. 


Figure  16.  Gate  sequencer. 


The  average-value  modeling  of  the  pulsed  power  weapon  is  based  on  the  6-pulse  load- 
commutate  converter  model  set  forth  previously  in  the  derivation  of  the  power  supply.  The 
significant  differences  are;  (1)  an  extension  to  the  12-pulse  case  is  not  necessary,  and  (2)  the 
representation  of  the  gate  sequence  control  is  different.  In  order  to  derive  an  average-value  model 
of  the  gate  sequence  controls  it  is  appropriate  to  begin  with  the  abc  model  of  the  filter.  In 
particular 

fl  -1  0^ 


dvj 


1 

0 


-1 


Uc  -  Vf 


(50) 


where  vyis  the  filtered  voltage  vector. 


Vf  =  [Vah  Vbc 


(51) 


Transforming  (51)  to  the  synchronous  reference  frame  yields 

1(3  Vs  ^ 

~  J 


'^‘if  _ 


dt 


-  Oie^df 


(52) 
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Table  7.  Pulsed  Power  Parameters. 

Name 

Description 

Value 

■ 

Charge  time 

0.8  s 

■i 

Catch-up  time 

0.1  s 

m 

Off  time 

0.05  s 

td 

Dead  time 

0.05  s 

c 

Pulse  Forming  capacitance 

4F 

L 

Pulse  forming  inductance 

3.21  mH 

r 

Resistance  at  pulse  forming  inductance 

2.42  mO. 

L. 

Transformer  leakage  inductance 

64.3  pH 

r, 

Transformer  resistance 

2.42  mO. 

Transformer  primary-to-secondary  turns  ratio 

5.00 

wm 

Current  gain 

0.3213  V/A 

Voltage  gain 

80A/V 

^lirn 

Current  limit 

6KA 

T 

Filter  time  constant 

4.594  ms 

Maximum  charging  voltage 

1000  V 

dv 


dt 


df  _ 


\(S 


V  - 


v,.„  +  —  V 


jp  ''df 


-  (Oy 


‘if 


(53) 


Together  (52)  and  (53)  are  used  to  represent  the  average-value  mode!  of  the  filter.  From  these 
equations  note  that  the  b  to  c  voltage  may  be  expressed 


Vhc=^V^-COS\ 


) 


(54) 


where 
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<t)j^^=mglc(v^+jvjf)  (55) 

from  which  it  can  be  seen  that  the  positive  zero  crossing  of  the  filtered  b  to  c  line-to-line  voltage 
occurs  at 

e  =--<!>.,  (56) 

Thus,  valve  3  is  fired  when  the  position  of  the  synchronous  reference  frame  is 


(57) 


By  definition,  valve  3  fires  at  an  angle  a  +  7tl3  after  the  peak  of  the  a-phase  primary  voltage. 
This  voltage  may  be  expressed  as 


%=V2v^cos(0,  +  <^>,,) 


(58) 


where 

=angle(v^-7v^) 

Therefore,  valve  3  fires  at 

/I  ^  A 

d,  =  a  +  --<t). 

Comparing  (60)  to  (57) 

~  ^ fpp  ~ 

The  average  value  model  of  the  pulsed  power  system  is  illustrated  in  Fig.  17. 


3 

Figure  17.  Average-value  model  of  pulsed  power  gate  sequencer. 


(59) 


(60) 


(61) 


27 


Harmonic  Filter 

In  order  to  reduce  the  harmonic  content  at  the  ac  bus,  a  wye-connected  series  RLC  filter  is 
used  (Fig.  18).  In  the  detailed  computer  model,  this  filter  is  represent  in  terms  of  abc  variables  as 


Figure  18.  Harmonic  filter. 


For  the  purpose  of  the  average-value  model,  it  is  convenient  to  represent  (62)  and  (63)  in  the 
synchronous  reference  frame 


The  filter  parameters  are  given  in  Table  8. 
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Hydrodynamic  Model 

The  hydrodynamic  model  used  in  the  representation  of  the  ship  propulsion  system  is 


_ T  \ 

dt  ~ 


where  Tg  is  the  induction  motor  torque  and  Tq  is  the  load  torque  represented  by 


(66) 

(67) 


The  assumed  value  of  parameters  are  listed  in  Table  9.  This  model  is  over  simplified  with  the 
inertia  selected  low  compared  to  the  actual  system.  A  more  detailed  model  can  be  readily 
substituted  for  the  model  presented  herein  if  the  actual  hydrodynamic  data  is  available. 


_ Table  9.  Hydrodynamic  Model _ _ 

7  =  50  ^Nms-  =  15.51  rad/s _ N=3 _ T^^  =  \.lMNm 


Component  Model  Interconnection 

Model  component  interconnection  is  always  an  issue  in  detailed  and  average-value 
modeling.  In  a  previous  SBIR,  PCKA  developed  interconnection  strategy  based  on  neglecting 
fast  transients  [5].  In  this  present  SBIR,  an  alternate  approach  of  retaining  those  transients  has 
been  explored.  This  is  accomplished  conveniently  by  representing  a  small  capacitance  at  the  ac 
bus.  Although  this  introduces  artificially  fast  time  constants,  these  fast  time  constants  do  not 
present  the  same  difficulties  which  are  introduced  in  the  case  of  detailed  simulations  since  the 
fast  time  constants  are  not  repeatedly  excited.  Therefore,  integration  algorithms  applicable  to 
stiff  systems  such  as  Gear's  Method,  are  effective.  In  terms  of  q-  and  d-variables,  the  capacitor 
is  represented  as 


C 


yv 


;=i 


(68) 


dt 


(69) 
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Equation  (68)  and  (69)  establish  the  ac  bus  voltage  as  state  variables.  The  capacitance  C  in  (68) 
and  (69)  was  selected  as  10  jiF. 
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COMPUTER  STUDIES 

Results  of  several  computer  studies  are  given  in  this  section  for  the  purposes  of 
demonstrating  the  facility  of  the  detailed  and  reduced-order  simulations  of  the  IPS  shown  in  Fig. 
1  and  to  compare  the  system  response  predicted  by  the  detailed  and  reduced-order  simulations. 
Four  studies  are  given.  With  the  exception  of  Study  2,  each  is  divided  into  two  parts.  In  part  "a" 
the  system  response  predicted  by  the  detailed  simulation  is  given;  in  part  "b"  the  predicted 
reduced-order  response  is  shown.  In  most  cases,  the  studies  are  self-explanatory;  therefore  a 
lengthy  discussion  is  not  given.  When  evaluating  the  studies,  it  is  convenient  to  refer  to  Fig.  1 
and  a  listing  of  the  variables  plotted  for  the  detailed  and  reduced-order  simulations.  For 
convenience  Fig.  1  and  the  list  of  variables  plotted  are  repeated  at  the  beginning  of  each  study. 
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Study  1 

In  this  study,  the  IPS  is  initially  operating  in  the  steady  state  with  the  ship  service  power 
supply  operating  at  full  load;  however,  the  pulsed  power  system  is  out-of-service.  At  t  =  15.5, 
the  torque  command  of  the  vector  controlled  15-phase  propulsion  induction  motor  is  increased 
from  zero  to  1.2  MN»m  (full  load).  The  dc  bus  stabilizer  is  incorporated  into  the  vector  control 
of  the  induction  motor. 

There  are  some  differences  between  the  system  responses  predicted  by  the  detailed  and  the 
reduced-order  simulations.  These  are  due  primarily  to  the  fact  that  PCKA  has  not  been  totally 
successful  in  developing  a  reduced-order  model  of  two  or  more  converters  operating  in  parallel 
from  the  same  bus.  Nevertheless,  the  accuracy  is  certainly  acceptable  for  this  application. 


Ship  Service  (DC) 
I  MW 


Figure  1.  System  studied. 
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Study  la 


Variable  Description  -  Detailed  Simulation 


ACSL  Name 

alphapwsp 

iabcimc(l) 

iabcppw(l) 

iabcpwsp(l) 

iabcsm(l) 

iaoiml(l) 

idcll 

idcllf 

iddrv 

idppw 

idpwsp 

idsm 

iqdrv 

iqppw 

iqpwsp 

iqsm 

irectpwsp 

irppw 

teiml 

tesm 

vabc(l) 

vaoiml(l) 

vcapppw 

vd 

vdcl  1 
vdcpwsp 
vfd 
vq 

wrmiml 


Description 

power  supply  rectifier  phase  delay  (rad) 

a-phase  current  into  induction  motor  power  converter  (A) 

a-phase  current  into  pulse  power  device  (A) 

a-phase  current  into  power  supply  (A) 

a-phase  synchronous  machine  current  (positive  in)  (A) 

a-phase  induction  motor  current  (A) 

dc  current  into  first  5  H-bridge  inverter  (not  filtered)  (A) 

low  pass  filtered  inverter  current  into  first  set  of  five  H-bridges 
(filtered  so  it  could  be  compared  to  average  value  model)  (A) 

d-axis  current  into  induction  motor  power  converter  (A) 

d-axis  current  into  pulse  power  device  (A) 

d-axis  current  into  power  supply  (A) 

d-axis  synchronous  machine  current  (positive  in)  (A) 

q-axis  current  into  induction  motor  power  converter  (A) 

q-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  power  supply  (A) 

q-axis  synchronous  machine  current  (positive  in)  (A) 

dc  link  current  in  power  supply  (positive  out  of  rectifier)  (A) 

inductor  current  in  pulse  power  device  (A) 

induction  motor  electromagnetic  torque  (Nm) 

synchronous  machine  torque  (positive  for  motor  operation)  (Nm) 

a-phase  4 1 60  V  bus  voltage  (V) 

a-phase  induction  motor  voltage  (V) 

capacitor  voltage  of  pulse  power  device  (A) 

d-axis  bus  voltage  (V) 

dc  voltage  at  input  of  first  5  H-bridges  (V) 

power  supply  dc  output  voltage  (V) 

synchronous  machine  field  voltage  (unreferred)  (V) 

q-axis  bus  voltage  (V) 

induction  motor  speed  -  mechanical  (rad/s) 


^ndDQA 
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Study  lb 


Variable  Description  -  Nonlinear  Reduced-Order  Model  Simulation 


ACSL  Name 

alphaspwsp 

idcll 

iddepspw(2) 

iqdeimc(l) 

iqdeimc(2) 

iqdepspw(l) 

iqdepwsp(  1 ) 

iqdepwsp(2) 

iqdesm(l) 

iqdesm(2) 

irectpwsp 

irppw 

teiml 

tesm 

vcapppw 

vdcl  1 

vdcpwsp 

vfd 

vqde( 1 ) 
vqde(2) 
wrmiml 


Description 

power  supply  rectifier  phase  delay  (rad) 

inverter  current  into  first  set  of  five  H-bridges  (A) 

d-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  induction  motor  power  converter  (A) 

d-axis  current  into  induction  motor  power  converter  (A) 

q-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  power  supply  (A) 

d-axis  current  into  power  supply  (A) 

q-axis  synchronous  machine  current  (positive  in)  (A) 

d-axis  synchronous  machine  current  (positive  in)  (A) 

dc  link  current  in  power  supply  (positive  out  of  rectifier)  (A) 

inductor  current  in  pulse  power  device  (A) 

induction  motor  electromagnetic  torque  (Nm) 

synchronous  machine  torque  (positive  for  motor  operation)  (Nm) 

capacitor  voltage  of  pulse  power  device  (A) 

dc  voltage  at  input  of  first  5  H-bridges  (V) 

power  supply  dc  output  voltage  (V) 

synchronous  machine  field  voltage  (unreferred)  (V) 

q-axis  bus  voltage  (V) 

d-axis  bus  voltage  (V) 

induction  motor  speed  -  mechanical  (rad/s) 


VDCll 

5000 


dsndi 
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Study  2 

In  this  study,  the  steady-state  waveforms  of  several  system  variables  are  shown.  The  ship 
service  power  supply  is  operating  at  full  load  with  the  pulsed  power  system  again  out-of-service. 
The  torque  command  for  the  propulsion  induction  motor  is  constant  at  full  load  (12  M  N»m). 


Ship  Service  (DC) 
IMW 


Figure  1 .  System  studied. 
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Variable  Description  -  Detailed  Simulation 


ACSL  Name 

alphapwsp 

iabcimc(l) 

iabcppw(l) 

iabcpwsp(l) 

iabcsm(l) 

iaoiml(l) 

idcll 

idcllf 

iddrv 

idppw 

idpwsp 

idsm 

iqdrv 

iqppw 

iqpwsp 

iqsm 

irectpwsp 

irppw 

teiml 

tesm 

vabc(l) 

vaoiml(l) 

vcapppw 

vd 

vdcl  1 
vdcpwsp 
vfd 
vq 

wrmiml 


Description 

power  supply  rectifier  phase  delay  (rad) 

a-phase  current  into  induction  motor  power  converter  (A) 

a-phase  current  into  pulse  power  device  (A) 

a-phase  current  into  power  supply  (A) 

a-phase  synchronous  machine  current  (positive  in)  (A) 

a-phase  induction  motor  current  (A) 

dc  current  into  first  5  H-bridge  inverter  (not  filtered)  (A) 

low  pass  filtered  inverter  current  into  first  set  of  five  H-bridges 
(filtered  so  it  could  be  compared  to  average  value  model)  (A) 

d-axis  current  into  induction  motor  power  converter  (A) 

d-axis  current  into  pulse  power  device  (A) 

d-axis  current  into  power  supply  (A) 

d-axis  synchronous  machine  current  (positive  in)  (A) 

q-axis  current  into  induction  motor  power  converter  (A) 

q-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  power  supply  (A) 

q-axis  synchronous  machine  current  (positive  in)  (A) 

dc  link  current  in  power  supply  (positive  out  of  rectifier)  (A) 

inductor  current  in  pulse  power  device  (A) 

induction  motor  electromagnetic  torque  (Nm) 

synchronous  machine  torque  (positive  for  motor  operation)  (Nm) 

a-phase  4160  V  bus  voltage  (V) 

a-phase  induction  motor  voltage  (V) 

capacitor  voltage  of  pulse  power  device  (A) 

d-axis  bus  voltage  (V) 

dc  voltage  at  input  of  first  5  H-bridges  (V) 

power  supply  dc  output  voltage  (V) 

synchronous  machine  field  voltage  (unreferred)  (V) 

q-axis  bus  voltage  (V) 

induction  motor  speed  -  mechanical  (rad/s) 


100000 


96/06/ 
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Study  3 


Initially  the  system  is  operating  in  the  steady  state  with  the  propulsion  system  (induction 
motor)  and  the  ship  service  power  supply  operating  at  full  load.  At  t  =  20.5  the  pulsed  power 
weapon  is  fired. 


Ship  Service  (DC) 
1  MW 


Figure  1 .  System  studied. 
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Study  3a 


Variable  Description  -  Detailed  Simulation 


ACSL  Name 

alphapwsp 

iabcimc(l) 

iabcppw(l) 

iabcpwsp(  1 ) 

iabcsm(l) 

iaoiml(l) 

idcll 

idcllf 

iddrv 

idppw 

idpwsp 

idsm 

iqdrv 

iqppw 

iqpwsp 

iqsm 

irectpwsp 

irppw 

teiml 

tesm 

vabc(l) 

vaoim  1(1) 

vcapppw 

vd 

vdcl  1 
vdcpwsp 
vfd 
vq 

wrmiml 


Description 

power  supply  rectifier  phase  delay  (rad) 

a-phase  current  into  induction  motor  power  converter  (A) 

a-phase  current  into  pulse  power  device  (A) 

a-phase  current  into  power  supply  (A) 

a-phase  synchronous  machine  current  (positive  in)  (A) 

a-phase  induction  motor  current  (A) 

dc  current  into  first  5  H-bridge  inverter  (not  filtered)  (A) 

low  pass  filtered  inverter  current  into  first  set  of  five  H-bridges 
(filtered  so  it  could  be  compared  to  average  value  model)  (A) 

d-axis  current  into  induction  motor  power  converter  (A) 

d-axis  current  into  pulse  power  device  (A) 

d-axis  current  into  power  supply  (A) 

d-axis  synchronous  machine  current  (positive  in)  (A) 

q-axis  current  into  induction  motor  power  converter  (A) 

q-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  power  supply  (A) 

q-axis  synchronous  machine  current  (positive  in)  (A) 

dc  link  current  in  power  supply  (positive  out  of  rectifier)  (A) 

inductor  current  in  pulse  power  device  (A) 

induction  motor  electromagnetic  torque  (Nm) 

synchronous  machine  torque  (positive  for  motor  operation)  (Nm) 

a-phase  4160  V  bus  voltage  (V) 

a-phase  induction  motor  voltage  (V) 

capacitor  voltage  of  pulse  power  device  (A) 

d-axis  bus  voltage  (V) 

dc  voltage  at  input  of  first  5  H-bridges  (V) 

power  supply  dc  output  voltage  (V) 

synchronous  machine  field  voltage  (unreferred)  (V) 

q-axis  bus  voltage  (V) 

induction  motor  speed  -  mechanical  (rad/s) 


VABC(l) 


61/90/96 


VDCPl 


96/06/ 
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Variable  Description 

ACSL  Name 

alphaspwsp 

idcll 

idclepspw(2) 

iqdeimc(l) 

iqdeiinc(2) 

iqdepspw(l) 

iqdepwsp(l) 

iqdepwsp(2) 

iqdesm(  1 ) 

iqdesm(2) 

irectpwsp 

irppw 

teiml 

tesm 

vcapppw 

vdcl  1 

vdcpwsp 

vfd 

vqde(  1 ) 
vqde(2) 
wrmiml 


Study  3b 

Nonlinear  Reduced-Order  Model  Simulation 
Description 

power  supply  rectifier  phase  delay  (rad) 

inverter  current  into  first  set  of  five  H-bridges  (A) 

d-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  induction  motor  power  converter  (A) 

d-axis  current  into  induction  motor  power  converter  (A) 

q-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  power  supply  (A) 

d-axis  current  into  power  supply  (A) 

q-axis  synchronous  machine  current  (positive  in)  (A) 

d-axis  synchronous  machine  current  (positive  in)  (A) 

dc  link  current  in  power  supply  (positive  out  of  rectifier)  (A) 

inductor  current  in  pulse  power  device  (A) 

induction  motor  electromagnetic  torque  (Nm) 

synchronous  machine  torque  (positive  for  motor  operation)  (Nm) 

capacitor  voltage  of  pulse  power  device  (A) 

dc  voltage  at  input  of  first  5  H-bridges  (V) 

power  supply  dc  output  voltage  (V) 

synchronous  machine  field  voltage  (unreferred)  (V) 

q-axis  bus  voltage  (V) 

d-axis  bus  voltage  (V) 

induction  motor  speed  -  mechanical  (rad/s) 


dsnd! 


IQDEPSPW(l) 


IRPPW 
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Study  4 

The  system  is  ,ni,ially  operating  with  the  propulsion  system  and  sh.p  service  power  supply 
operating  at  full  load  with  the  stabilising  control  incorporated  into  the  vector  control  of  the 
propulsion  induction  motor.  At  r  =  24.5  the  pulsed  power  weapon  is  fired  and  at  the  same  time 
the  stabilizing  control  is  deactivated  whereupon  the  propulsion  tnduction  motor  ,s  controlled  by 
the  "standard-  vector  control.  During  the  operation  of  the  weapon,  the  instab, hty  due  to  the 
standard  vector  control  is  noted  in  the  wavefomrs.  Compare  TESM  in  Study  3  to  TESM  .n 
Smdy  4  However,  once  the  weapon  has  discharged  the  instability  is  more  or  less  present  only  .n 
the  dc  variables  as  Illustrated  by  the  plots  of  VDCll  and  IDCl  IF.  Note  that  with  the  stab.hzer, 

before  the  finng,  the  vanables  are  stable;  however,  after  the  finng.  when  the  stabrlizer  .s  ou.-of- 

.  V  -.vchahnitv  Tn  Studv  3  wherein  the  stabilizer  was  in  service, 

service,  the  variables  demonstrate  instability.  In  btudy  , 

the  variables  return  to  the  pre-firing  stable  waveforms. 

In  the  case  of  the  reduced-order  simulation,  the  plots  stop  at  approximately  r  =  25.4.  This  ,s 

due  to  the  fact  that  the  reduced-order  model  detects  that  a  commutation  fa, lure  occurred  and 

subsequently  stopped  computing  since  thereafter  the  reduced-order  model  is  no  longer  vahd. 


Figure  1 .  System  studied. 
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Study  4a 


Variable  Description  -  Detailed  Simulation 


ACSL  Name 

alphapwsp 

iabcimc(l) 

iabcppw(l) 

iabcpwsp(l) 

iabcsm(l) 

iaoiml(l) 

idcll 

idcl If 

iddrv 

idppw 

idpwsp 

idsm 

iqdrv 

iqppw 

iqpwsp 

iqsm 

irectpwsp 

irppw 

teiml 

tesm 

vabc(  1 ) 

vaoiml(l) 

vcapppw 

vd 

vdcl  1 
vdcpwsp 
vfd 
vq 

wrmiml 


Description 

power  supply  rectifier  phase  delay  (rad) 

a-phase  current  into  induction  motor  power  converter  (A) 

a-phase  current  into  pulse  power  device  (A) 

a-phase  current  into  power  supply  (A) 

a-phase  synchronous  machine  current  (positive  in)  (A) 

a-phase  induction  motor  current  (A) 

dc  current  into  first  5  H-bridge  inverter  (not  filtered)  (A) 

low  pass  filtered  inverter  current  into  first  set  of  five  H-bridges 
(filtered  so  it  could  be  compared  to  average  value  model)  (A) 

d-axis  current  into  induction  motor  power  converter  (A) 

d-axis  current  into  pulse  power  device  (A) 

d-axis  current  into  power  supply  (A) 

d-axis  synchronous  machine  current  (positive  in)  (A) 

q-axis  current  into  induction  motor  power  converter  (A) 

q-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  power  supply  (A) 

q-axis  synchronous  machine  current  (positive  in)  (A) 

dc  link  current  in  power  supply  (positive  out  of  rectifier)  (A) 

inductor  current  in  pulse  power  device  (A) 

induction  motor  electromagnetic  torque  (Nm) 

synchronous  machine  torque  (positive  for  motor  operation)  (Nm) 

a-phase  4160  V  bus  voltage  (V) 

a-phase  induction  motor  voltage  (V) 

capacitor  voltage  of  pulse  power  device  (A) 

d-axis  bus  voltage  (V) 

dc  voltage  at  input  of  first  5  H-bridges  (V) 

power  supply  dc  output  voltage  (V) 

synchronous  machine  field  voltage  (unreferred)  (V) 

q-axis  bus  voltage  (V) 

induction  motor  speed  -  mechanical  (rad/s) 


96/06/ 
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Study  4b 


Variable  Description  -  Nonlinear  Reduced-Order  Model  Simulation 


ACSL  Name 

alphaspwsp 

idcll 

iddepspw(2) 

iqdeimc(l) 

iqdeimc(2) 

iqdepspw(l) 

iqdepwsp(  1 ) 

iqdepwsp(2) 

iqdesm(l) 

iqdesm(2) 

irectpwsp 

irppw 

teiml 

tesm 

vcapppw 

vdcl  1 

vdcpwsp 

vfd 

vqde(l) 

vqde(2) 

wrmiml 


Description 

power  supply  rectifier  phase  delay  (rad) 

inverter  current  into  first  set  of  five  H-bridges  (A) 

d-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  induction  motor  power  converter  (A) 

d-axis  current  into  induction  motor  power  converter  (A) 

q-axis  current  into  pulse  power  device  (A) 

q-axis  current  into  power  supply  (A) 

d-axis  current  into  power  supply  (A) 

q-axis  synchronous  machine  current  (positive  in)  (A) 

d-axis  synchronous  machine  current  (positive  in)  (A) 

dc  link  current  in  power  supply  (positive  out  of  rectifier)  (A) 

inductor  current  in  pulse  power  device  (A) 

induction  motor  electromagnetic  torque  (Nm) 

synchronous  machine  torque  (positive  for  motor  operation)  (Nm) 

capacitor  voltage  of  pulse  power  device  (A) 

dc  voltage  at  input  of  first  5  H-bridges  (V) 

power  supply  dc  output  voltage  (V) 

synchronous  machine  field  voltage  (unreterred)  (V) 

q-axis  bus  voltage  (V) 

d-axis  bus  voltage  (V) 

induction  motor  speed  -  mechanical  (rad/s) 


dsnd! 
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APPENDICES 

The  following  appendices  would  normally  be  included  as  sections  of  the  report;  however, 
each  appendix  is  written  in  manuscript  format  for  publication  purposes.  It  was  decided  to 
include  these  manuscripts  as  part  of  the  report  and  not  to  spend  the  time  to  rewrite  this  material 
in  the  format  of  the  report. 

Appendix  A;  DC  Link  Stabilized  Field  Oriented  Control  of  Electric  Propulsion  Systems 
Appendix  B:  Modeling  Consideration  in  a  15-Phase  Induction  Motor  Drive  System 
Appendix  C:  Control  of  a  15-Phase  Induction  Motor  Drive  System 


Appendix  A 

DC  Link  Stabilized  Field  Oriented  Control  of  Electric  Propulsion  Systems 


Abstract  -  In  the  process  of  designing  the  controller  for  the 
induction  motor  drive  included  in  this  report,  it  was  discovered 
that  the  system  would  exhibit  negative  impedance  instability. 
Therefore,  a  simple  though  nonlinear  control  algorithm  was 
created  to  mitigate  this  Instability.  As  a  fint  step  in  developing 
this  control  for  the  IPS  system,  the  control  was  developed  for  a 
3-phase  induction  motor  based  drive  system.  In  particular,  the 
3-phase  induction  motor  based  version  of  the  control  was 
verified  both  through  the  use  of  computer  simulation  and  by 
construction  of  a  3.7  kW  prototype  system  nsing  the  laboratory 
facilities  of  the  University  of  Missouri  -  RoUa.  The  design 
and  performance  of  the  3-phase  version  of  the  control  are 
set  forth  herein. 

L  Introduction 

Electric  propulsion  systems  consisting  of  a  turbine 
driven  synchronous  machine  feeding  an  induction  motor 
drive  through  a  rectifier  -  dc  link  -  inverter  fiequency 
changer  are  important  in  a  wide  variety  of  applications 
including  locomotives,  hybrid  electric  vehicles,  submarines, 
and  ship  propulsion  systems,  including  the  next  generation 
warship,  the  Surfece  Combatant  21.  The  advantage  of  such 
a  system  over  a  Tn»*chi»nical  transmission  is  that  the  tuibine 
speed  becomes  fully  decoupled  from  the  load  speed,  allowing 
the  tuibine  speed  to  be  optimized  with  regard  to  fuel 
efficiency.  In  addition,  the  elimination  of  the  mechanical 
linkage  between  the  turbine  or  other  prime  mover  and  the 
mechanical  load  (drive  train  or  propeller)  allows  a  greater 
degree  of  architectural  freedom  in  the  locomotive  /  vehicle  i 
ship  layout  In  many  of  these  systems,  the  inverter  is  used  to 
tightly  regulate  the  motor  current  waveforms,  which  has  the 
advantage  of  making  the  inverter  /  motor  drive  extremely 
robust  with  regard  to  overcurrents.  However,  at  the  same 
timi-!  such  regulation  has  the  disadvantage  that  it  makes  the 
motor  drive  appear  to  have  a  negative  impedance,  since  if 
the  inverter  voltage  is  reduced  the  dc  link  cunent  will 
increase  so  as  to  maintain  constant  power  (since  the  motor 
currents  will  remain  undisturbed).  This  negative  impedance 
can  result  in  loss  of  dynamic  stability  of  the  propulsion 
system.  In  order  to  avoid  instability,  a  classical  method  is  to 
increase  the  dc  link  capacitance.  However,  in  large  drive 
systems,  such  as  those  used  in  naval  applications,  the 
capacitance  required  can  become  costly  in  terms  of  capitaL 
space,  and  weight.  In  this  paper,  a  nonlinear  dc  link 
stabilized  field  oriented  control  is  demonstrated  which  is 
shown  have  performance  characteristics  similar  to  the 
classical  field  oriented  control,  but  mingaies  the  negative 
impedance  instability  problem. 


DL  System  Overview 

Fig.  1  illustrates  the  type  of  electric  propulsion  system 
considered  herein.  The  power  source  of  this  system  is  a 
djif^i  '“nginff  or  turbine  which  serves  as  a  prime  mover  for 
the  3-phase  synchronous  machine  (SM)-  The  3-phase  output 
of  the  machine  is  rectified  using  an  uncontrolled  rectifier. 
The  rectifier  output  voltage  is  denoted  Vr.  An  LC  circuit 
serves  as  a  filter,  and  the  output  of  this  fiilter  is  denoted  v*,  • 

A  voltage  regulator  /  exciter  adjusts  the  field  voltage  of  the 
synchronous  machine  in  such  a  way  that  the  source  bus 
voltage  Vdc3  is  equal  to  the  commanded  bus  voltage  The 
source  bus  is  connected  via  a  tie  line  to  the  load  bus,  the 
voltage  at  which  is  denoted  v^ci  *  The  load  bus  consists  of  a 
capacitive  filter  (which  includes  both  electrolytic  and 
polypropylene  capacitance)  as  well  as  a  3-phase  fully 
controlled  inverter,  which  is  in  turn  supplies  an  induction 
motor.  The  induction  motor  drives  the  mechanical  load, 
which  is  rotating  at  a  speed  Based  upon  the 

mw-haniral  rotor  spccd,  and  the  desired  electromagnetic 
torque  (which  is  determined  by  the  controller 

governing  the  mechanical  system),  the  induction  motor 
controls  specify  the  on/off  status  of  each  of  the  inverter 
semiconductors  in  such  a  way  that  the  desired  torque  is 
obtained.  Although  this  system  is  quite  robust  with  regard  to 
overcurrents,  and  simple  to  design  from  the  viewpoint  that 
the  controller  governing  the  mechanical  system  is  decoupled 
from  the  control  of  the  electrical  system  (since  the  torque  can 
be  controller  nearly  instantaneously),  such  systems  are  prone 
to  be  subjea  to  a  limit  cycle  behavior  in  the  dc  bus  voltage 
known  as  negative  impedance  instability. 

m.  Cause  Of  Negative  Impedance  Instability 

In  order  to  gain  insight  into  the  cause  of  negative 
impedance  instability,  it  is  appropriate  to  set  forth  a  highly 
simplified  model  of  the  system  depicted  in  Fig.  1.  This 
model  will  focus  on  the  dc  bus  dynamics  and  need  only  be 
valid  in  the  tens  to  hundreds  of  hertz  band.  An  appropriate 
model  of  the  synchronous  machine  /  rectifier  and  LC  filter  is 
illustrated  in  Fig.  2.  Therein,  the  dynamics  of  both  the 
prime  mover  and  voltage  regulator  are  neglected  since  these 
arc  subject  to  long  time  constants  imposed  by  the  prime 
mover  inertia  and  synchronous  machine  field,  respecttvely. 
The  dc  voltage  behind  inductance  and  resistance 
synchronous  machine  model  is  based  on  the  work  set  forth  in 
[11,  which  has  been  shown  to  have  excellent  bandwidth.  In 
Fig.  2, 
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Figure  1 .  Electric  propulsion  system. 


£  =  (1) 

where  oorjm  is  the  electrical  rotor  speed  of  the  synchronous 
marhinp  and  and  A,"  are  the  q-  and  d-axis  subtransient 
flux  linifages,  a  is  the  firing  angle  relative  to  the 
subtransient  back  cmt  ic(P)  and  L,(p)  are  the  conunutating 
and  transient  commutating  inductance's  defined  by 

Z,c(P)  =  +(Ij  -i")sin(2p  +f)  (2) 

and 

L,(P)  =  L"<i  +L"  +  iL"  -  L")sin(2p  -  f)  (3) 


where  it  is  assumed  that  the  instantaneous  power  is  equal  to 
the  instantaneous  power  command  P*  defined  as 

P'  =  (4) 

In  (4),  7^  is  an  instantaneous  torque  command  which  is  the 
input  to  the  field  oriented  induction  motor  control. 
Typically,  the  instantaneous  torque  command  T,  is  set 
equal  to  the  desired  torque  T,,du-  However,  the  control 
aignrithm  proposed  in  this  paper  will  possess  an  alternate 
relationship. 

Upon  n(>glt»rting  the  tie  line,  which  is  generally  short  in 
an  electrical  sense  for  the  frequency  range  of  interest,  the 
component  models  illustrated  in  Fig.  2.  and  Fig.  3.  are 
combined  as  in  Fig.  4.  Therein, 


where  P  is  the  firing  angle  relative  to  rotor  position  and  L" 
and  L"  are  the  synchronous  machine  q-  and  d-axis 
subtransient  inductance’s,  is  the  synchronous  machine 
stator  resistance,  Lf  and  /y  are  the  inductance  and  resistance 
of  the  LC  filter  inductor,  and  Cf  is  the  capacitance  of  the  LC 
filter  capacitor. 

A  highly  simplified  model  of  the  capacitive  filter  and 
induction  motor  is  illustrated  in  Fig.  3.  Therein,  the  drive  is 
modeled  as  a  dependent  current  source  equal  to  the 
instantaneous  power  P  divided  by  dc  link  voltage  vja  and 
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Synchronous  Machine/Rectifier  Filter 


L,  =  L,m+Lf 

(5) 

fZ(c(P)ti5rsm  "T  2/* sjm 

(6) 

C,  =  Cf+Ci 

(7) 

’„  =  ^^^£cos  (a) 

(8) 

In  (5-8),  a  and  p  can  be  found  by  steady-state  analysis  of  the 
load  commutated  converter  rectifier  system  as  set  forth  in 
[2].  However,  aggnming  the  rectifier  is  uncontrolled  the 
r^ioilatinn  may  be  made  much  simpler  by  neglecting 
subtransient  saliency  whereupon 


Figure  2.  Simplified  model  of  synchronous 
machine/rectifier  and  filter. 


Figure  3.  Simplified  capacitive  filter, 

induction  motor  drive  model. 
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a  =  0 


(9) 


and  the  commutating  and  transient  commutating 
inductance's  are  no  longer  a  function  of  P .  Furthermore, 
due  to  the  action  of  the  voltage  regulator,  the  subtransient 
inductance's  will  have  a  value  such  that 


V«  =  v^. 


(10) 


so  that  no  subtransient  information  is  actually  needed.  It 
should  be  cautioned  that  this  model  is  intended  for 
explanation  purposes  and  for  guidance  in  designing  control 
algorithms,  not  for  high-fidelity  simulation  or  for  the  testing 
of  control  algorithms. 

In  order  to  utilize  this  equivalent  circuit  to  predict 
negative  impedance  stability,  note  that  linearizing  the  input 
current  with  respect  to  input  voltage  about  an  operating  point 
wherein  Vdd  is  equal  to  Vda  yields 
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from  which  it  is  apparent  that  the  inverter  presents  an  input 
impedance  of 
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which  has  a  characteristic  equation  of 


M  J  1 


P’Re^ 

'  c,lJ 


=  0 


From  (14)  it  follows  that  necessary 

conditions  for  stability 

are  that 

(15) 

and  that 

R. ,  P’ 

(16) 

I.  C.v’A. 

Physically  (15)  is  a  steady  state  limit  on  the  amount  of 
power  which  can  be  transferred  and  is  normally  satisfied  in 
practice.  In  the  case  of  (16),  it  can  be  seen  that  a  given 
system  will  eventually  go  unstable  as  the  power  command  is 
increased.  To  avoid  instability,  one  method  is  to  simply 
increase  C,  until  the  system  is  stable  to  the  maximum 
possible  power  command.  However,  this  has  the 
disadvantage  in  that  for  very  large  drive  systems  the  physical 
space,  weight,  and  cost  of  this  capacitance  become  issues, 
espedaiiy  considering  the  fact  that  the  applications  are 
largely  mobile  in  nature. 

IV,  Lipoc  Stabilizing  Field  -  Oriented  Control 

In  this  sectioiL  an  algorithm  which  improves  the 
damping  of  the  dc  link  by  eliminated  the  negative  impedance 
ffffprt  over  a  prescribed  bandwidth  is  set  forth.  This 
algorithm  is  based  upon  the  fea  that  torque  control  in  a  field 
oriented  drive  is  nearly  instantaneous.  As  mentioned 
previously,  typically  the  instantaneous  torque  command  T,  is 
set  equal  to  the  desired  torque  determined  by  the 

control  algorithm  governing  the  mechanical  system. 
However,  herein  it  is  proposed  to  determine  the 
instantaneous  torque  command  as 


Te,des 


As  can  be  seeiL  in  a  small  signal  sense  the  inverter  appears 
to  possess  a  negative  impedance  which,  from  an  intuitive 
point  of  view,  would  seem  to  be  destabilizing. 

In  order  to  verify  this  conclusion,  replacing  the 
dependent  source  representing  the  inverter  with  its  small 
signal  equivalent  impedance,  yields  the  following  small 
signal  model  of  the  equivalent  circuit 


where  v*,  is  the  filtered  dc  inverter  voltage,  i.e., 

dvjci  Vdd  -Vdci 
dt  ^ 


(16) 


(17) 


(13) 


and  the  parameters  and  n  and  x  are  considered  to  be 
constants  herein  but  could  also  be  made  to  be  a  function  of 
operation  point 

The  advantage  of  this  simple  though  nonlinear  control 
algorithm  is  that  it  is  extremely  straightforward  to 
implement  yet  highly  effective  in  mitigating  negative 
impedance  instabilities.  In  order  to  illustrate  the  effect  of  the 
algorithm  on  the  system  note  that  using  the  control  law  input 
power  into  the  inverter  is  given  by 


(14) 


(18) 


where 


P des  ~  Pe.des^ 


(19) 


From  (18)  the  input  current  may  be  expressed 
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Linearizing  (20)  about  the  desired  operating  pomt(vrf«  = 

)  yields 

=  (21) 


n-  IPdMs 


If  T  so  great  as  to  not  interact  with  the  dc  link  dynamics  and 
n  is  selected  to  be  unity  then  the  input  impedance  presented 
by  the  inverter  is  infinite  over  the  frequency  range  over 
which  negative  impedance  instabilities  occur,  thus  avoiding 
this  type  of  instability. 

Although  the  explanation  of  the  previous  paragraph 
illustrates  the  basic  philosophy  of  the  control,  the 
possibilities  of  the  control  are  much  richer  than  is  indicated 
therein.  In  particular,  by  suitable  selection  of  t  and  n  a 
large  variety  of  behaviors  can  be  obtained.  In  order  to  see 
this  it  is  helpful  to  first  set  forth  the  nonlinear  differential 
equation  governing  the  dc  link  dynamics  in  the  presence  of 
the  new  control.  In  particular. 
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Linearization  of  (22)  yields 
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In  order  to  illustrate  the  effects  of  varying  nand  t, 
consider  the  case  of  a  system  in  which  =  400  V, 
i?,  =  4.58Q.  I.  =  13.9  iiiH,  and  C/=51.4  nF.  These 
parameters  correspond  to  a  test  system  which  was  used  for 
laboratory  demonstration.  Fig.  5  illustrates  the  root  ioci  of 
the  characteristic  equation  as  t  is  varied  from  0.1  ms  to  1  s 
(24)  for  «=i.3,5,  and  7  (It  should  be  noted  that  n  does  not 
have  to  be  an  integer).  As  can  be  seen,  in  each  case  the  root 
locus  contains  an  unstable  complex  pole  (denoted  A  and  A*) 
for  values  of  t  which  becomes  stable  as  x  is  increased. 
For  all  n  shown  in  Fig.  5  the  real  part  of  the  eigenvalues 


becomes  more  negative  as  x  in  increased.  In  addition, 
initially  the  complex  part  also  decreases.  In  the  case  of  /» =  5 
eventually  the  complex  pair  becomes  real  (point  B)  and  then 
one  of  these  real  roots  meets  the  root  corresponding  the  filter 
at  point  C,  at  which  this  pair  of  eigenvalues  becomes 
complex.  In  the  case  of  ;i  =  7  the  two  complex  poles 
eventually  become  real  at  point  D;  after  which  the  pair 
moves  away  from  each  other  on  the  real  axis. 

Fig.  6  illustrates  the  damping  of  the  complex  pole  pair 
as  n  and  x  are  varied.  Note  that  for  each  value  of  n  there  is  a 
value  of  xwhich  maximizes  the  damping.  It  is  also  apparent 
that  generally  speaking,  as  n  is  increased  the  damping  can 
be  increased. 

Since  the  control  law  is  nonlinear,  it  will  of  course  be  a 
function  of  operating  poinc  and  so  it  is  important  to 
investigate  the  performance  of  the  control  as  the  operating 
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Figure  7.  Root  locus  versus  P  and  n. 

point  (primarily  through  power  command)  is  varied.  Fig.  7 
illustrates  the  root  locus  of  the  system  as  power  command  is 
varied  with  (n  =  1,  t  =  4  ms),  (n  =  3,  t  =  2.4  ms),  (n  =  5,  t  = 
2.7  ms),  and  (n  =  7,  X  =  3.1  ms).  In  each  case  x  was  selected 
so  as  to  maximize  the  damping  factor  of  the  complex  pole 
pair.  As  can  be  see  from  Fig.  7,  at  low  power  the  complex 
eigenvalue  of  the  system  is  at  point  (A,A*)  irregardless  of  n . 
In  the  case  of  (n=l),  the  location  of  the  roots  is  power  level 
independent  However,  in  the  case  of  /i=3  the  complex  pair 
moves  to  (B,B*).  Finally,  in  the  case  of  n=5,  and  n=l,  the 
complex  eigenvalues  become  real  at  point  C  and  D, 
respectively.  This  illustrates  an  important  feature  of  the 
proposed  control  law  which  is  that  although  with  the 
standard  control  (rt=0)  the  system  becomes  less  stable  as 
power  level  increases  with  the  proposed  control  law  the 
system  actually  becomes  mote  stable  as  power  level 
increases,  with  the  e,xception  of  (n=l)  in  which  case  the 
eigenvalues  associated  with  the  dc  link  become  largely 
operating  point  invariant. 

V.  Implementation 

Before  setting  forth  the  implantation  of  the  proposed 
controller,  it  is  appropriate  to  first  consider  a  standard  field 
oriented  control  such  as  the  rotor  flux  indirea  field  oriented 
control  illustrated  in  Fig.  8  (note  that  the  control  proposed  in 
this  paper,  is,  however,  independent  of  whether  or  not  the 
field  oriented  control  is  direa  or  indirect).  Therein,  an 
instantaneous  torque  command  is  the  input  to  the 
controller.  This  torque  command  is  equal  to  the  torque 
desired  by  the  controller  governing  the  mechanical 
dynamics,  As  can  be  seen,  based  on  the  torque 

command  H  and  desired  d-axis  rotor  flux  level  Xj’,  the 
desired  q-  and  d-  axis  stator  currents,  and  /*,  are 
determined.  This  calculation  is  a  function  of  the  induction 


motor  rotor  magnetizing  inductance  ,  the  induction  motor 
rotor  inductance  (rotor  leakage  plus  magnetizing)  I^,,  the 
rotor  resistance  ,  and  the  number  of  poles,  Np^m  ■  Based 
on  the  q-  and  d-  axis  stator  currents  the  electrical  radian  slip 
frequency,  coj,,„,  is  determined,  which  is  then  added  to  the 
electrical  rotor  speed  cDr.im  in  order  to  determine  the  electrical 
speed  of  the  synchronous  reference  frame  which  is 

integrated  in  order  to  determine  the  position  of  the 
synchronous  reference  firame  In  addition  to  the 

algorithm  illustrated  in  Fig.  8,  especially  in  large  drives,  the 
field  oriented  control  will  often  include  on  line  parameter 
iHpntifiratifin  algorithm  to  Compensate  for  variations  of  the 
rotor  time  constant 

Once  the  q-  and  d-aids  current  commands  and  the 
position  of  the  synchronous  reference  finme  are  established. 
thp^  currents  may  be  synthesized  in  a  variety  of  ways. 
Herein,  the  q-  and  d-axis  current  command  was  transformed 
back  into  a  abc  variable  current  command  which  is  an  input 
to  a  hysteresis  type  current  control. 

Incorporating  the  link  stabilizing  control  into  the  field 
oriented  control  is  quite  straightforward.  In  particular,  the 
only  rfiffhrpTira;  in  the  control  is  the  instantaneous  torque 
command  is  generated  using  (16)  rather  than  being  set  equal 
to  the  desired  torque,  as  is  illustrated  in  Fig.  9. 

VL  Experimental  setup 

In  order  to  illustrate  the  experimental  set  up  a  system 
such  as  the  one  depiaed  in  Fig.  1  was  constructed  at  a  low 
(3.7  kW)  power  level.  The  prime  mover  was  a  dynamometer 
in  speed  control  mode.  The  parameters  of  the  3.7  kW 
synchronous  machine  are  listed  in  Table  1.  ThereirL  all 
rotor  parameters  have  been  referred  to  the  stator  by  the 
appropriate  turns  ratio.  A  solid  state  e.xciter  voltage  regulator 
was  used  to  control  the  dc  link  voltage;  a  block  diagram  of 


Figure  8.  Rotor  Flux  Oriented  Indirect  Field  Oriented 
Control. 


edes" 


^dci’ 


^dcumin 


^dcLmax 

_ z:  _ 

J 

1 

D 

T  S+\ 

Figure  9.  DC  link  stabilizing  control. 

this  control  appears  in  Fig.  10.  The  LC  filter,  transmission 
line,  and  inverter  capacitance  parameters  are  listed  in  Table 
2.  Finally^  a  3.7  kW  induction  motor  was  used  as  a  load. 
Induction  motor  parameters  are  listed  in  Table  3.  The 
mHiirrinn  motor  contTol  was  implemented  based  on  the 
indirect  rotor  flux  strategy,  and  the  current  command  was 
synthesized  using  hysteresis  current  control  with  a  hysteresis 
band  of  0.95  A.  The  link  stabiUzing  controller  parameters 
were  n  =  1,  -c  =  4  ms,  =  200  V,  and  =  600  V. 


Table  1.  Synchronous  machine  parameters 

r  =  382  mW 

r.  sm 

i,,=0.83  mH 

P  =  4 

L  =  13.5  mH 

ma 

L^^=39.7mH 

n  ^,.^=0.0271 

r^,  =  31.8n 

122  mW 

=  0.923  n 

Z,,^  =  2.54  mH 

'•m/=  40.47  0 

^Ikdl  “ 

r,.,,=  1.3in 

=  3.68  mH 

VIL  Validation 

The  performance  of  the  link  stabilizing  field  oriented 
control  was  validated  using  both  a  detailed  (as  opposed  to 
average  value  /  reduced  order  )  computer  simulation  and  in 
hardware.  For  the  purposes  of  computer  simulation,  the 
synchronous  machine  and  induction  motor  models  used  were 
those  set  forth  in  [3].  In  the  case  of  the  salient  pole 
synchronous  machine,  magnetic  saturation  was  represented 


^,=  1 


Table  2.  Passive  component  parameters. 

L^=  4.69  mH 

r^=  1.54  n 

C^=  10.1  |iF 

r,.  =  273  m£2 

!me _ 

C,  =  41.3  pF 

Table  3.  Induction  motor  parameters. 

r  =  400  mfi 

5  _  _ 

L,^  =  5.73  mH 

L  =64.3mH 

r/  =  227  mfi 

Z,,/  =  4.94tnH 

?  =  4 

in  the  d-axis.  The  simulation  included  the  switching  of  each 
power  electronics  device.  Semiconductor  conduction  losses 
were  inciudei  however  sv/itching  losses  were  neglected. 

Figure  11  illustrates  the  simulated  performance  of  the 
system  as  the  desired  torque  is  changed  from  2  to  19  Nm 
over  a  period  of  100  ms.  Variables  depicted  included  the 
commanded  a-phase  current  the  actual  a-phase  current 
/ox,  the  dc  inverter  voltage  v^dy  and  the  electromagnetic 
torque  Although  the  actual  torque  closely  tracks  the 
desired  torque,  it  can  be  seen  that  as  the  power  command 
increases  the  dc  bus  voltage  becomes  unstable,  stressing  both 
the  semi-conductors  and  the  capacitors.  In  a  typical  system, 
such  behavior  could  easily  result  in  the  semiconductor  and/or 
capacitor  failure.  The  experimental  system  has  been 
constructed  so  as  to  be  able  to  survive  the  overvoltages.  Fig. 
12  depicts  the  same  study  as  measured  in  the  laboratory.  In 
Fig.  12,  the  instantaneous  electromagnetic  torque  is  not 
shown  because  suitable  instrumentation  was  not  available. 
As  can  be  seen  there  is  a  reasonable  correspondence  between 
Fig.  11  and  Fig.  12  with  the  exception  that  the  actual  system 
appears  to  be  less  stable  than  is  predicted  by  the  simulation. 
This  is  because  the  power  requirements  of  the  actual  drive 
system  are  greater  than  the  simulated  system  because  of 
switching  losses  (It  should  be  noted  that  the  magmtude  of 
the  voltage  swing  increases  very  rapidly  with  power  level). 
In  addition,  once  a  system  becomes  unstable  it  tends  to  be 
very  sensitive  to  parameter  variations. 

Fig.  13  depicts  the  performance  of  the  same  system 
with  the  link-  stabilizing  field  oriented  control  as  calculated 
using  the  computer  simulation.  As  can  be  seen,  according  to 
the  simulation  the  torque  still  closely  tracks  the  commanded 
torque,  however  in  this  case  there  is  no  evidence  of 
instability.  Fig.  14  illustrates  the  system  performance  as 
measured  in  the  laboratory.  As  predicted,  the  dc  bus  voltage 
is  well  behaved  and  the  dc  link  bus  voltage  is  stable. 

One  concern  which  may  arise  is  a  possible  reduction  in 
torque  bandwidth  since  a  drop  in  inverter  voltage  will  result 
in  a  transient  dip  in  torque.  The  detailed  computer 
simulation  was  used  to  investigate  this  effect  since  the 
primary  variable  of  interest  was  the  electromagnetic  torque. 
Fig.  15  depicts  the  predicted  change  of  performance  of  the 
standard  field  oriented  control  as  the  torque  command  is 
stepped  from  2  to  19  Nm.  As  can  be  seen,  the 
electromagnetic  torque  reaches  the  commanded  value  in 


Figure  10.  Voltage  regulator/exciter. 


Figure  i  1 .  Simulated  performance  of  standard  fieW 
oriented  control  during  ramp  increase  in 
desired  torque. 


Figure  13.  Simulated  performance  of  link  stabilized 
field  oriented  control  during  ramp  increase 
in  desired  torque. 


Figure  12.  Measured  performance  of  standard  field 
oriented  control  during  ramp  increase  in 
desired  torque. 

approximately  5  ms.  The  torque  response  is  not 
instantaneous  due  to  the  fact  that  a  step  change  in  current 
cannot  be  achieved  in  practice  and  because  the  dip  in  link 
voltage  causes  a  temporary  loss  of  current  tracking  in  the 
hysteresis  current  control.  Fig.  16  depicts  the  response  of  the 
link  stabilized  field  oriented  control.  In  this  case,  the 
electromagnetic  torque  reaches  the  commanded  value  in  the 
order  of  8  ms.  Although  the  link  stabilized  control  is 
somewhat  slower  than  the  standard  field  oriented  control, 
this  sLieht  reduction  in  bandwidth  is  not  a  significant 
disadvantage  in  view  of  the  improved  dc  bus  voltage.  This  is 
particuiariy  true  in  the  fact  that  most  propulsion  systems 


Figure  14.  Measured  performance  ot  link  stabilized 
field  oriented  control  during  ramp 
increase  in  desired  torque. 

have  mechanical  inertia  such  that  in  either  case  the  torque 
response  may  be  considered  to  be  instantaneous. 

vnL  Conclusions 

A  straightforward  but  nonlinear  control  algorithm  has 
been  set  forth  which  can  be  used  to  mitigaie  negative 
impedance  instabilities.  The  effectiveness  of  the  control  has 
been  demonstrated  both  through  the  use  of  computer 
simulation  and  in  the  laboratory,  and  found  to  be  highly 
effective  in  eliminating  negative  impedance  instability.  The 
primary  disadvantage  of  the  control  is  a  slightly  slower 
torque  response.  In  addition  to  being  applicable  to  induction 


Figure  15.  Simulated  performance  of  standard  field 
oriented  control  during  step  change  in 
desired  torque. 


Figure  16.  Simulated  performance  of  link  stabilized 
field  oriented  control  during  step  change 
in  desired  torque. 

motor  drives,  the  control  algorithm  could  also  be  used  with 
other  types  of  machines  in  which  rapid  torque  control  is 
possible,  most  notably  permanent  magnet  synchronous 
machines. 
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Appendix  B 

Modeling  Consideration  in  a  15-Phase  Induction  Motor  Drive  System 


L  Introduction 

This  analytical  work  set  forth  herein  is  the  basis  of  both 
the  detailed  and  average  value  models  used  in  this  SBIR 
Phase  I  effort.  The  work  has  been  documented  in  journal 
format  so  that  it  may  be  published  in  the  technical  literature 
at  a  future  date. 


BL  System  Description 

Fig.  1  illustrates  the  inverter  topology  of  the  15-phase 
induction  motor  drive  system  which  is  the  subject  of  this 
analysis.  As  can  be  seen,  the  inverter  is  comprised  of  three 
dc  rails,  each  of  which  is  conneaed  to  five  H-bridge 
inveners.  These  inverters  arc  designated  H,  where  the  i‘th 
H-bridge  supplies  the  /*th  phase  of  the  machine.  The  phases 
of  the  machine  are  in  turn  labeled  such  that  the  fundamental 
component  of  the  i'th  phase  is  displaced  (i-l)24^  (electrical) 
from  the  1st  phase. 

The  voltage  vector  v^cwill  be  used  to  designate  the 
individual  rail  voltages,  i.e., 

Vrfc  ^  Vrfci  Vffcj  ]  (^) 


and  similarly  the  current  veaor  irfcwili  be  used  to  designate 
the  dc  link  currents 


idc 


-[ 


^clc\ 


^dci 


(2) 


regulates  the  current  through  a  voltage  command  rather  than 
by  a  direct  method  such  as  hysteresis  current  control.  The 
focus  of  this  work  is  on  the  basic  model  of  the  machine  and 
inverter.  Therefore,  the  discussion  will  concentrate  on  the 
representation  of  only  the  innermost  control  loop,  which  is 
based  on  the  sine-triangle  modulation  strategy  in  which  a 
sinusoidal  duty  cycle  reference  is  compared  to  a  triangle 
modulation  waveform,  as  illustrated  in  Fig.  3.  The 
switching  signal  for  the  i'th  phase  is  then  synthesized  as 


I  1  ifsi*>  di 
[  -ItfSj  ^di 


(4) 


The  zero  state  is  not  being  used  in  order  to  reduce  common 
mode  link  currents. 


m.  Machine  Variable  Induction  Machine  Model 


The  machine  variable  induction  motor  consists  of  the 
stator  and  rotor  voltage  equations,  the  flux  linkage 
equations,  and  the  torque  equation.  The  stator  and  rotor 
voltage  equations  of  a  symmetrical  N-phase  induction 
machine  may  be  expressed  as 


and 


(5) 


The  configuration  of  each  H-bridge  is  illustrated  in  Fig. 
2.  Therein  each  semiconductor  actually  represents  a  number 
of  devices  in  series,  as  necessary  to  achieve  the  required 
voltage  blocking  capability,  and  £>i,  through  Z)4,  designate 
the  on/off  status  of  each  of  the  active  semiconductors  of  the 
;'th  bridge,  and  Vj,,  i„,  and  /'dc*,  denote  the  i'th  phase 
voltage,  the  i'th  phase  current  (defined  positive  into  the 
machine),  and  the  dc  link  current  into  the  i'th  H-bridge, 
respectiveiy.  The  switching  state  of  the  i'th  H-bridge  is 
defined  as 


1  D\,  and £>4,  on 

0  (£>i,  on  and  Dt,  on)  or  {Di,  on  and  Da,  on)  (3) 

-1  D2,  and  £>3,  on 


and  s  will  denote  the  vector  of  switching  states. 

The  drive  system  under  consideration  will  be  controlled 
using  the  field-oriented  principal.  The  current  command 
will  be  achieved  using  a  current  control  algorithm  that 


‘del 


‘del 


‘dci 


Figure  I .  Drive  topology. 


(11) 


Figure  2.  H-bridge. 


respectively,  wherein 


“[/u  fis  ...  /ns  ] 


(6) 

O) 


and 


where  Lu  is  the  stator  leakage  inductance,  I//  is  an  N  by  N 
identity  matrix,  L  ms  is  the  stator  magnetizing  inductance,  and 
C«  is  a  N  by  N  matrix  whose  elements  are  given  by 


Css,j=cos{(j-i)f)  (12) 


The  stator-to-rotor  seif  inductance  matrix  Lsr 
expressed 

=  Cjr  Ssr  ] 

In  (13),  the  elements  of  Cjrand  SirOre  defined  by 


may  be 


(13) 


C^=cos(9r-(/-l)f)  (14) 

and 

S,r,  =Sm(0r-(/-l)f)  (15) 


respectively,  Lsr  denotes  the  maximum  mutual  inductance 
between  any  stator  and  any  effective  rotor  winding,  and  9  r  is 
the  electrical  rotor  position,  defined  as  PH  times  the 
mechanical  rotor  position,  Qrmy  where  P  is  the  number  of 
poles. 

In  (10),  the  rotor  self  inductance  is  of  the  form 


f  r  —  [  far  fbr  J  (8) 

where  f  can  be  a  voltage  v,  flux  linkage  A,  or  current  i.  As 
can  be  seen,  the  distributed  rotor  circuits  (which  are 
physically  implemented  with  a  squirrel  cage  type 
construction)  are  represented  on  a  2-phase  basis  regardless  of 
the  number  of  stator  phases. 

Upon  neglecting  the  effects  of  magnetic  saturation  and 
space  harmonics,  the  stator  and  rotor  flux  linkage  equation 
may  be  expressed 

\s  =  Lsjij  +  Ljrir  (9) 

Ar  =  LS.i,+L,rir  (10) 

In  (9),  the  stator  self-inductance  matrix  is  given  by 


Lrr-  ^  {Lir+L  mr)l2  ( 16) 

where  hr  and  Lmr^c  the  rotor  leakage  and  magnetizing 
inductances,  respectively. 

From  the  flux  linkage  equations,  the  expression  for 
electromagnetic  torque  is  readily  established  using  coenergy 
techniques.  In  particuiar.  in  terms  of  machine  variables 

Te  =  Ssr  Ssr  ]  Jr  (17) 


It  should  be  noted  that  Lmr,  and  Z.,^arc  not 
independent  because  they  share  the  same  magnetizing  path. 
In  particular. 


L  -  — 

L.ms  -  „ 

(18) 

L 

(19) 

T  risfir  • 

(20) 

where  n,znd  rir  are  the  effective  number  of  turns  of  the 
stator  and  rotor  windings,  respectively,  and  R„  denotes  the 

reluctance  of  the  magnetizing  path. 


Figure  3.  Sine-triangle  modulation  signals  d-  and 


IV.  Referred  Machine  Variable  Model 

The  next  step  in  deriving  the  qd  model  for  the  induction 
motor  is  to  refer  the  rotor  variables  to  the  stator  by  the 
appropriate  turns  ratio.  In  particular,  it  is  convenient  to 
define  the  referred  rotor  voltage  and  flux  linkage  as 

(21) 


The  referred  rotor  current  is  defined  as 


i'  -  ±.iLL: 
•r  ~ 


where  is  an  integer  which  is  dependent  upon  the 
construction  the  qd  transformation  (this  will  be  discussed  in 
a  following  section). 

In  terms  of  referred  variables,  the  voltage  equations 


become 


v,  =  r..,  +  — 

y/  +  ^ 


where  the  referred  rotor  resistance  given  by 


The  referred  flux  linkage  equations  are  given  by 


inverse  are  also  constant  matrices  which  further  promotes 
computationally  efficiency.  In  the  case  of  average-value 
modeling,  selection  of  a  synchronous  reference  frame  is  most 
convenient  as  such  a  formulation  will  result  in  state 
variables  which  are  constant  in  the  steady-state.  Such  a 
description  is  readily  linearized.  Furthermore,  if  the 
dynamics  of  the  system  in  which  the  drive  is  place  are  stiff,  a 
model  in  which  the  state-variables  are  constant  in  the 
steady-state  may  be  very  efficiently  solved  using  Gear's 
method. 

However,  beyond  the  choice  of  reference  frame  for 
detailed  and  average-value  model,  in  the  case  of  drive 
systems  with  a  large  number  of  phases  the  construction  of 
the  qd  transformation-  which  is  not  unique,  becomes  an 
issue.  For  example,  for  the  15-phase  system  considered 
herein,  the  transformation  could  be  made  with  a  single  q- 
and  d-axis  and  13  zero  sequence  variables,  two  q-  and  d- 
axis  with  11  zero  sequence  variables,  etc.  As  it  turns  ouL 
two  parricuiarly  interestmg  choices  are  the  use  of  five  sets  of 
q-  and  d-axis  with  five  zero  sequence  variables,  and  three 
sets  of  q-  and  d-axis  with  nine  zero  sequence  variables, 
which  correspond  to  the  partitioning  the  machine  into  five 
3-pfaase  sets  and  three  5-phasc  sets,  respectively. 

The  construction  of  a  qd  transformation  based  on  five 
3-phase  sets  is  straightforward  to  construct  because  on  a 
3-phase  set  basis  it  is  identical  to  the  classical  3-phase  qd 
transformation.  In  this  case  the  transformation  is 
constructed  by  dividing  the  fifteen  phases  into  five  3-phase 
sets  and  arranging  the  elements  of  the  conventional  3 -phase 
qd  transformation  accordingly.  In  particular,  the  stator  and 
rotor  transformations  may  be  expressed  as 


Ai  —  Li5jij  +  Lsrir 

(27) 

t\r  =  L/i,  +  L(rir 

(28) 

In  (25) 

=  |£«[  C,r  -S,r  ] 

(29) 

and 

Lrr  = 

(30) 

where 

(31) 

V.  QD  Transformations 

The  transformation  of  the  referred  machine  variable 
model  to  a  qd  model  has  several  advantages.  First  of  all,  in 
the  case  of  detailed  computer  simulation,  such  a 
transformation  eliminates  rotor  position  dependence  in  the 
flux  linkage  equations  so  that  the  inverse  inductance  matrix 
only  has  to  be  computed  once  in  the  event  that  flux  linkages 
are  used  as  state  variables.  By  choosing  the  stationary 
reference  frame,  the  reference  fiame  transformation  and  its 


^qds  “  Kffy 

and 

(32) 

(33) 

respectively.  In  (32) 

iqds  ”  I  ^qds\  ^qdsr  ^qds-i  ^qdsA 

(34) 

where 

fgdsi  =  L  ] 

(35) 

In  (33) 

^qdr  —  fqr  fdr  J 

The  stator  transfonnation  matrix  is  given  by 

(36) 

!  K,„  ••• 

K,  =  |  : 

i  K,,, 

(37) 

where  ail  the  elements  of  K„j  are  zero  except  for  the  for  the 
following: 

(K,„),,  =  |cos(e  -(5(/--  1) +»■  -  1)§)  (38) 

(K,,,)z,  =  |sin(9  -(5(/-  - 1)  +  «  - 1)^)  (39) 

(K,„)3,  =  |  (40) 

where©  is  the  position  of  the  arbitrary  reference  frame.  The 
rotor  transformation  matrix  is  given  by 


Kr  = 


cos(9-0r)  sin(9-0r) 
sin(9  -  9  r)  -cos(9  -  0  r) 


(41) 


Transformation  of  the  stator  and  referred  rotor  voltage 
equations  (with  ot=3)  yields 

dkq,, 

V qsi  “  ^ qsf  (OKdsi  ^ 


dt 

+  CO A,<7j,  H 


VOj,  +* 


dt 


and 


dka 


Vqr  =  rriqr  +  {(0-(i)  r)^dr  + 
Vdr  =  r'ridr  -  (cO  -  ©  r)^,r  + 


(42) 

(43) 

(44) 

(45) 

(46) 


Finally^  the  electromagnetic  torque  may  be  expressed 


Tt  —  2  O^dmiqsi  ^qmidai) 

L  L  y-\ 


(54) 


The  qd  equivalent  circuit  corresponding  to  (42-52)  is 
illustrated  in  Fig.  4.  As  can  be  seen  the  equivalent  circuit  is 
comprised  of  a  q-axis  circuit  d-axis  circuit  and  five 
identical  and  non-coupled  zero-sequence  circuits. 


i  ,  r,  L,, 


4* 


Vi 


^qs3 


+ 

^qs3 


^qs4 


4* 

^qx4 


^qsS 


^r^dsS 
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/r  '  r 


respectively,  where  /  e  1, 5 .  Transformation  of  the  referred  ^<*2 


dsl  ^ 


yields 

“  d^ls^qsi  ^qrn 
~  Lls^ds{  4  X-t/zTi 

(47) 

(48) 

'^dxi 

-7 - 

‘ds4 

and 

Kqr  =  Lj^iqr 
^dr  ^  ^ir^dr  ^ dm 

(49) 

(50) 

^Js4 

'^ds5 

where 

^qm  ~  Ziml  2  Iqsf  4/qr  j 

(51) 

"^ds5 
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(52) 
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Figure  4.  Five  3-phase  set  QD  equivalent  circuit  of 
1 5-ohase  induction  machine. 


The  3 -phase  set  based  equivalent  circuit  has  several 
attractive  characteristics.  First,  the  transformation  is 
essentially  identical  to  the  standard  3-phase  transformation, 
at  least  on  a  per  set  basis.  Second,  the  zero  sequence  circuits 
are  identical  and  decoupled.  In  addition,  the  equivalent 
circuit  is  very  suggestive  from  a  control  point  of  view  to 
hanrilP!  the  situation  in  which  one  H-bridge  is  disabled.  For 
example,  suppose  that  the  drive  is  used  in  the  context  of  a 
field-oriented  control  system  and  a  fault  is  detected  on  one 
H-btidge  corresponding  to  one  phase  of  the  /’th  3-phase 
stator  set.  Then  the  i'th  q-  and  d-axis  current  command 
would  be  set  equal  to  zero,  and  the  other  two  H-bridges 
with  the  3 -phase  set  would  also  be  disabled.  Next 
the  total  d-axis  stator  current  command  would  be  maintained 
by  increasing  the  d-axis  current  command  in  the  other  for 
sets,  so  as  to  maintain  constant  rotor  flux.  Depending  upon 
the  torque  command  the  q-axis  current  commands  in  the 
other  phases  could  be  increased  so  as  to  maintain  constant 
torque,  unless  the  current  command  has  reached  its  limit  in 
which  case  the  torque  output  would  be  forced  to  decrease. 
Even  so,  this  strategf  has  advantages  in  that  machine 
operation  would  be  balanced  which  would  avoid  unexpeaed 
thermal  losses  and  torque  harmonics. 

Another  useful  choice  of  reference  frames  is  based  on  a 
5-pbase  set  based  approach  with  three  q-  and  d-axis  and  nine 
zero  sequence  variables.  This  choice  is  pamculariy 
convenient  because  such  a  construction  associates  groups  ail 
the  phases  conneaed  to  a  given  dc  link  into  a  set  Before 
considering  the  full  15  by  15  transformation,  it  is  appropriate 
to  set  forth  the  transformation  of  a  5-phase  set  to  q-  and 
d-axis  variables.  Such  a  transformation  may  be  expressed  as 


where 

and 


fsj  =  [  /u  ■■■  fu  ] 


fqd5s  —  ^  yijj  /ds  Xos  fas  Jbs  ] 


(55) 

(56) 

(57) 


In  order  to  derive  the  construction  of  Ksj ,  it  is  convenient  to 
express  it  in  terms  of  rows  as 


qi 


K5,= 


qs 


(38) 


where  q,  denotes  the  i'th  row  of  Ksj.  The  first  and  second 
rows  of  Ksi  formulate  the  q-  and  d-axis.  From  geometrical 
arguments,  the  first  two  rows  may  be  expressed  as 


qi  = 

cos(0)  cos(9 


f)  cos(0-y)  cos(0-f)  cos(0-f)  ] 

(59) 


sin(0)  sin(0-x)  sin(0-Y)  sin(0-y)  sin(0-Y)  ] 

(60) 

It  is  convenient  to  define  the  third  row  of  the  transformation 
as 

q3  =  [HHT]  (61) 


With  this  definition,  /ojis  the  instantaneous  average  phase 
current  in  the  5-phase  set.  Furthermore  (61)  is  clearly 
orthogonal  to  both  (59)  and  (60).  Definition  of  the  other 
sequence  variables,  /ouand  /pj,  is  somewhat  more  difficult 
because  it  is  not  clear  what  the  physical  interpretation,  if 
any,  should  be.  Herein,  the  row  vectors  q4and  qswill  be 
obtained  using  the  Gram-Schmidt  algorithm,  in  which  a 
non-ortho  gonai  set  of  basis  vectors  is  transformed  to  a 
orthogonal  one.  The  non-orthogonal  set  is  assumed  to  consist 
of  qi,q2,q3.a4,  and  as,  where,  as  previously  noted,  qi 
through  q3  are  already  mutually  orthogonal  and 


a4  =  [10000] 
a5  =  [oi0  0o] 


(62) 

(63) 


which  vieids 

:i/2 


^2  -I-Zcosy  -I-Zcosy  -1-2 


cosy  -1-2 COS 


t] 

(64) 


and 
q5  = 


5  J 150 +  100  cos  f 

[0  2-2cosy  -1-4cosy  -1-4cosy  -3-2cosy  ] 

^  (65) 

Thus  the  5 -phase  based  transformation  is  given  by 
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2 
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-l-2cos^ 

2j2  ^ 

-I--4cos‘y 

2j2 

-1-4  cos  ^ 

2V2 

-3-2 cos ^ 

Jl5(KI00cca^ 

JlSMOOcosy 
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The  extension  of  the  5-phase  transformation  to  the 
15-phase  case  begins  with  the  definition  of  the  full  15-phase 
stator  transformation.  In  particular. 


foots  —  Kjfj 


where 

and 


^qds  ~  ^  fqrfu  J 

fqdsi  ~  1^  Jqsi  fdsi  fos!  fasi  f^i  J 


fth  coiumn(Ks,)  = 


fcos(e-(/-i)f^) 

fsin(e-(/-l)i) 

[K5sl4.  J-M 
[K5*]5^ 


set  forth  in  (67).  The  stator  and  rotor  voltage  equations 
become 

Vqsf  =  ^s^qsi  -rCOAociy,  i - r--  (73) 


d^dsf 


^ dsi  ^ dsf  ^^qsf 


'^Ojf  ~  ^s^Ost  ' 


VrtSf  —  ^slcLSt  ' 


V  pjf  —  r  si  + ' 


“  fri qr  r^^dr 

Vdr  =  r'ridr  -  (CO  -  (0  r)Xqr  + 


It  should  be  noted  that  the  same  symbolism  is  being  used  as 
in  the  previous  15 -phase  transformation.  However,  this 
should  not  arise  because  the  two  transformations  are  never 
used  at  the  same  time.  In  this  case,  the  transformation  Kjis 
formed  so  as  to  apply  K5,  to  three  independent  5-phase  sets 
(this  will  require  adjusting  the  phase  of  the  reference  frame 
position  dependent  terms).  This  yields 

'k,,  ‘ 

K,=  K,,  (70) 

where  all  elements  of  K„  are  zero  except  those  in  the  i'th. 
3-^/'th,  6+/’th,  9+/’th,  and  12+/*th  columns.  In  particular,  the 
j’th  column  of  K„  is  given  by 


where  /  e  1,2, 3 .  The  flux  stator  linkage  equations  may  be 
expressed  as 

^^Isiqsi  ^  '^qm  (^0) 

^dsf  ^^Isidsi  ‘^dm  (^I) 

AOji  ^Lisiosi  (^2) 

~  ^IsictSi 

y^^si-Lisi^si  (84) 

and  the  rotor  flux  linkage  equations  as 


/^qr  —  ^Isiqr  ^qm 

(85) 

A/^r  ~  ^isidr  "^dm 

(86) 

^qm  ~  m  1  2  /'  qsf  i  qr  j 

7  3  ^ 

(87) 

m  S  idsf  i  dr  J 

(88) 

L,m  ~  T-^ms 

(89) 

where 


Finally^  electromagnetic  torque  may  be  expressed 

Te  ~  iqsi  ^  idst  )  (^0) 

The  three  5-pbase  set  based  qd  equivalent  circuit  of  the 
induction  machine  is  illustrated  in  Fig.  5. 


y  e  /,  3  +  /,  6  +  /,  9  4-  /,  12 + /  (72)  VL  Detailed  Inverter  Model  And  Drive  Simulation 

The  rotor  transformation  is  identical  to  that  given  by  (33)  detailed  simulation  of  the  inverter  drive  makes  use  of 

and  (41).  ^  machine  variable  inverter  model  and  the  3-phase  set  based 

The  next  step  is  to  actually  transform  the  referred  model  in  the  stationary  reference  frame.  The  state 
variable  machine  model  (with  m-5)  using  the  transformation  variables  of  the  electrical  system  will  be  the  stator  and  rotor 

flux  linkages.  The  simulation  structure  is  illustrated  in  Fig. 
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Figure  5.  Three  5-phase  set  QD  equivalent  circuit  of 
1 5-phase  induction  machine. 

6.  The  algorithm  used  to  calculate  the  derivative  of  the  flux 
linkages  based  on  the  input  and  state  variables  is  as  follows. 
First,  given  the  flux  linkages,  the  qd  currents  are  found  by 
solving  (47-52).  From  the  qd  current  vector  the  machine 
currents  are  found  using  the  inverse  transformation  (which  is 
constant  if  the  stationary  reference  frame  is  selected,  and 
sparse).  Next  the  switching  state  veaor  is  found  in 
accordance  with  (3).  Based  on  the  switching  vector  s  and  the 
machine  current  i,the  voltage  vector  can  be  found.  In 
particular. 
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Figure  6.  Detailed  simulation  structure. 

Vdcj  - 2v(  if s,  =  1  and 
Vdcj  +  2va  if  Sj  =  1  and  /„  <  0 
^  -Vie,  +  2vt  if  s,  =  -1  and  /„  <  0 
-Vdc)  -  2vd  if  s*  =  -1  and  is,  >  0 


(91) 


where  j  corresponds  to  the  dc  link  of  the  i’th  phase.  The  dc 
current  into  the  i'th  H-bridge  is  given  by 


is,  if  Sf  =  1 
-is,  ifj,  =  -l 


(92) 


Finally,  the  dc  link  current  may  be  expressed 

5 

i  dc,  ^  i  dchsit~\y*.j  (93) 

i=i 


where  j  e  1,2,3.  Based  on  the  machine  variable  voltage 
vector  determined  by  (91)  the  qd  vanables  whereupon 
(42-46)  may  be  used  to  calculate  the  time  derivative  of  the 
flux  linkages. 

VTL  Average- Value  Lnverter  Model  And  Drive 
Simulation 

The  average-value  model  presented  herein  is  based  upon 
the  5-phase  set  based  qd  model  using  a  synchronous 
reference  frame.  The  first  step  necessary  in  order  to 
formulate  the  complete  model  is  to  determine  the  input 
voltage  veaor  Vqds  (the  line  indicates  fast  average  value) 
based  upon  the  dc  voltage  veaor  Vdc  and  the  qd  duty  c\’cle 
veaor  ddq  ■  (obtained  by  transforming  d  to  qd  variables).  For 
the  normal  situation  in  which  the  sine  triangle  algorithm  is 
not  overmodulated  dqd  is  such  that 

\di\  <  1  /  e  1...15  (94) 

it  can  be  readily  shown  (since  the  fast  average  of  the  output 
voltage  is  directly  proportional  to  the  input  voltage)  that 

^qds,  -  VdCjddqj  j  S  1...3  (95) 


In  the  overmodulated  case,  which  is  to  say  that 


ua>i  (96) 

for  some  /,  then  assuming  the  instantaneous  overmodulation 
is  due  to  the  q-  and  d-axis  components  rather  than  the  zero 
sequence  components  (which  is  normally  the  case  since  the 
zero  sequence  components  are  normally  zero)  then  using 
Fourier  analysis  to  extract  the  fundamental  component  it  can 
be  shown  that 

Vqd,  =  |)  J  ^  ^  (9'^ 

where 


The  structure  of  the  average  value  is  illustrated  in  Fig.  7. 
Therein,  all  qd  variables  and  electromagnetic  torque  should 
be  interpreted  on  a  fast  average  basis,  and  the  qd  machine 
model  is  the  three  5-phase  set  based  model  in  the 
svTichronous  reference  frame. 


Figure  7.  Average-value  simulation  structure.. 


Once  the  qd  input  voltage  has  been  established,  the  qd 
machine  model  may  be  used  to  determine  the  qd  current 
vector  and  calculate  the  derivative  of  flux  linkages.  It  can 
also  be  used  to  calculate  the  electromagnetic  torque,  as  is 
necessary  to  compute  the  mechanical  dynamics.  In  order  to 
interfece  the  drive  simulation  to  the  rest  of  the  electncal 
system  (the  dc  supplies)  one  further  calculation  is  necessaiy-. 
that  is  the  calculation  of  the  dc  link  currents.  This  is  based 
on  power  balance. 

In  the  case  of  the  five-phase  transformation  (66),  it  can 
be  shown  that  the  power  into  a  five-phase  set  is  given  by 

P  =  xiVqsiqs  +  Vdsids  +  2{VOsios  +  Vasias  +  Vp,/  p,))  (1001 

Therefore  the  power  into  the  j'th  dc  link  is  given  by 

Pj  =  \{Vqsjiqs,+Vds,ids,  +2{Vosposj  +Vasj‘asj +V9s)i^))) 

(101) 

The  power  into  the  j'th  link  may  also  be  expressed  as 

Pj  =  Vdcjidcj  (^02) 

Equating  the  two  results  and  incorporating  (95)  yields 

idcj  —  ^{dqsiqs  -rddsids  +2(doijl(li  •i'dasias  +  d^si^))  (lOl'i 


or,  in  the  case  of  overmodulated  operation 
'dci  —  ^(dqsiqs  +  dds^ds  +  2(dosios  +  das^as  + 


(104) 


Appendix  C 

Control  of  a  15-Phase  Induction  Motor  Drive  System 


L  Introduction 

The  work  set  forth  herein  has  been  conducted  for  Ae 
SBIR  Phase  I.  In  particular,  in  order  to  include  the  IPS  drive 
system  load  a  substantial  effort  had  to  be  placed  into  the 
control  design.  The  woric  has  been  formatted  in  journal 
format  so  that  it  may  be  published  in  the  technical  literature 
at  a  future  date. 


EL  System  Description 

The  application  of  the  drive  system  under  consideration 
is  a  Naval  propulsion  system  illustrated  in  Fig.  1.  Therein,  a 
turbine  is  the  prime  mover  for  a  3-phase,  60-Hz,  20  MW, 
4160  V  1-1  rms  synchronous  machine,  whose  output  voltage 
is  conuolled  bv  an  IEEE  Type  2  excitation  system.  The 
principal  load  on  the  bus  is  a  19  MW  induction  motor  drive 
system,  although  the  synchronous  machine  also  provides 
ship  service  power.  Parameters  of  the  synchronous  machine 
and  exciter  are  listed  in  Table  1. 


Fig.  2  illustrates  the  inverter  topology  of  the  15-phase 
induction  motor  drive  system  which  is  the  subject  of  this 


Figure  1.  System  Architecture. 


analysis.  As  can  be  seen  from  Fig.  1,  the  motor  power 
converter  is  divided  into  three  parallel  rectifier  -  dc  link 
inverter  paths.  In  each  path,  a  6-pulse  uncontrolled  rectifier 
output  is  filtered  with  an  LC  circuit  (I^=  3  mH  and  C^=  1 
mF),  the  output  of  which  is  coimected  to  five  H-bridge 
inverters,  wherein  Hi  designates  the  /'th  H-bridge  which 
controls  the  applied  to  the  /'th  phase  of  the  machine.  The 
phases  of  the  machine  are  labeled  such  that  the  fundamental 
component  of  the  voltage  appUed  to  the  /’th  phase  is 
displaced  (i-l)24^  (electrical)  from  the  first  phase. 

The  voltage  vector  v<fcwill  be  used  to  designate  the 
individual  rail  voltages,  i.e., 

Vrfc  =  [  ^dcx  ^dcz  ^dCi  J 

An  identical  convention  will  be  used  to  identify  the  dc 
current  vector  idt,  the  rectifier  current  vector  ir,  and  the 
rectifier  voltage  vector  Vr . 

The  configuration  of  each  H-bridge  is  illustrated  in  Fig. 

3.  Therein,  each  semiconductor  shown  represents  an 
effective  device  which  is  physically  constructed  usmg  a 
series  connection  of  several  devices  in  order  to  achieve  die 
necessary  voltage  rating.  The  on/off  status  of  each  effective 
switching  device  in  the  i*th  bridge  is  designated  Di,  through 
,  and  v,, ,  ,  and  idch,  denote  the  /’th  phase  voltage,  /’th 

phase  current  (defined  positive  into  the  machine),  and  the  dc 
link  current  into  the  /’th  H-bridge,  respectively.  The 
switching  state  of  the  i’th  H-bridge  will  be  denoted 

1  Z)i,  and£)4,  on 

Si  =  ^  0  (Dx^on  and  Da,  on)  or  (Da,  on  and  D4,  on)  (2) 
-1  D2,  and  D3,  on 

snd  s  will  denote  the  vector  of  switching  states.  The  goal  of 
the  controller  set  forth  herein  will  be  to  specify  s  as  a 
function  of  time  in  such  a  way  that  the  electromagnetic 
torque  T,  closely  tracks  the  desired  torque  H  • 

m.  Induction  Motor  Model 
The  control  of  the  induction  motor  is  formulated  on  the 
three  5-phase  set  based  qd  transformation  as  set  forth  in 
Appendix  B.  In  particular,  the  stator  transformation  may  be 
expressed  as 

fqdst  ~  (^) 

where 


In  (6) 

9„=0,-(/-l)^  (7) 

where  0^  is  the  position  of  the  synchronous  reference  frame. 
The  rotor  transformation  is  given  by 

fqdr  =K.fr  (8) 

where 


Figure  3.  H-bridge. 


(9)  The  three  5-phase  set  based  qd  equivalent  circuit  of  the 
induction  machine  is  illustrated  in  Fig.  4,  and  parameters  are 
listed  in  Table  2. 


In  (11)  0r  denotes  the  electrical  rotor  position.  In  terms  of 
qd  variables,  the  induction  motor  stator  and  rotor  voltage 
equations  may  be  expressed  as 


and 


^qst  “  ^ qsi  ® 

-  . 

Vfisf  ~ 

VOs,  =  Tj/Oy,  + 

Volt/  =  ^5^0?/  “• 

dk^Si 

vp,,  =rwp.,+-^ 


^qr  =  fiiqr  +  (Pe  - 

=  ^r^dr  ^  (®  e  “  ®  r)kqr 


respectively.  Since  the  machine  is  of  squirrel  cage  design, 
the  q-  and  d-axis  rotor  voltages  are  both  zero.  The  stator 
flux  linkage  equations  may  be  expressed  as 


^qst  ~Lls^qSi 

(19) 

'^dsi  ~^ls^dsi 

(20) 

Ao5f  “  Lisiosi 

(21) 

k<xsj  ~  Ligiasi 

(22) 

A-pjf  ^Lisi^si 

(23) 

and  the  rotor  flux  linkage  equations  as 

A-gr  —  ^Is^qr  ^qm 

(24) 

=  Llsifir  kdm 

(25) 

where 

kqm  —  Zfnl  2  Iqsj  +  t qr  1 

(26) 

kdm  ”  Lm\  2  idsi  "^^dr) 

\f=l  ' 

(27) 

Finally,  electromagnetic  torque  may  be  expressed 

Te  —  ^qsi  ^  kqm^  /<&,  j  (28) 


*•  r  CO,.  A.-,/ 


ielr5 


Figure  4.  Three  5-phase  set  QD  equivalent  circuit  of 
1 5-phase  induction  machine. 


IV.  Controller  Archttecture 
The  first  decision  to  be  made  in  the  design  of  the 
controller  was  the  overall  control  architecture.  Some  of  the 
desired  properties  of  the  controller  included  that  (i)  it  must 
be  stable,  (//)  it  is  field  oriented  based  so  that  the  drive  will 
act  as  a  torque  transducer,  (iii)  the  currents  be  tightly 
controlled  in  order  to  protect  the  inverter  semiconductors, 
and  (/v)  the  switching  flrequency  should  be  a  constant.  In 
regard  to  (/),  most  systems  are  inherently  stable.  However, 
because  of  requirement  (iii),  that  is  the  currents  are  to  be 
tightly  controlled,  this  drive  system  exhibits  a  negative 
impedance  characteristic  which  makes  the  drive  inherently 
unstable.  Requirement  (if)  is  put  in  place  because  is 
simplifies  the  system  design  in  that  the  with  a  field-oriented 
control  the  response  of  the  torque  response  of  the  drive  is 
much  faster  than  the  response  of  the  mechanical  dynamics, 
so  that  in  designing  the  mecbamcal  system  controls  the 
electrical  dynamics  may  be  neglected.  In  addition,  the  fast 
action  of  the  field  oriented  control  will  be  useful  in 
stabilizing  the  dc  link,  which  as  previously  mentioned  is 
unstable  if  the  motor  currents  are  tightly  regulated,  which  is 


in  turn  a  requirement  designed  to  protect  the  inverter 
semiconductors  (by  tightly  controlling  the  currents  the 
semiconductors  can  be  operated  close  to  there  maximum 
current  limit,  which  is  important  in  large  drive  systems). 
Finally  requirement  (iv)  is  put  in  place  so  that  switching 
losses  can  be  controlled  independently  of  operating  point. 

These  four  constraints  suggest  the  controller  architecture 
illustrated  in  Fig.  5.  The  structure  of  the  controls  closely 
parallels  that  of  the  system  in  that  there  are  three  conttol 
rhanneU,  one  for  each  dc  link.  The  primary  input  commands 
are  operation  status,  a  torque  command  and  a  rotor  flux 
command  The  first  stage  of  the  control  is  a  paralleling 
control  which  divides  the  torque  command  into  three  torque 
commands  corresponding  to  the  amount  of  electromagnetic 
torque  is  due  to  the  q-axis  current  associated  with  each  link, 
7^, ,  7^2 .  and  7^3 ,  respectively.  The  paralleling  controls  also 
divide  the  d-axis  rotor  flux  into  the  amount  of  that  flux 
which  is  to  be  contributed  by  the  5-phase  set  being  supplied 
by  that  link  These  commands  are  designated  X,  j,., ,  ,  and 

X  .  For  normal  operation  this  process  is  very  simple,  as 
illustrated  in  Fig.  6.  In  the  case  that  say,  one  dc  link 


Fieiire  5.  1 5-Pha.se  induction  motor  link  stabilized  indirect  field  oriented  controller. 


Figure  6.  Paralleling  controls. 


becomes  non-operationai  the  torque  and  flux  command  is 
divided  between  the  two-remaining  dc  links.  The  next  block 
in  the  controls  is  a  torque  command  limiter,  as  illustrated  in 
Fig.  7.  Note  that  this  limit  not  only  bounds  the  torque 
command  but  also  the  rate  of  change  of  the  torque  command, 
if  desired.  The  outputs  of  the  torque  limiters  are  the  link 

torque  commands  denoted  ,  2^2 »  ^3  • 

Based  on  the  link  flux  linkage  commands  and  the  link 
torque  commands,  the  desired  q-  and  d-axis  currents  for  each 
link  /‘,,and  /*,  are  established  using  the  link  stabilizing 
control  illustrated  in  Fig.  8.  The  function  of  this  control  is  to 
modify  the  qd  current  commands  such  that  the  negative 
impedance  instability  is  avoided.  This  is  accomplished  by 
making  the  inverter  appear  to  be  an  open  circuit  at 
flequencies  much  greater  than  1/x.  Details  are  set  forth  in 
Appendix  A.  The  outputs  of  the  link  stabilizing  control  are 
then  passed  through  a  current  limiter,  as  illustrated  in  Fig.  9. 
Note  that  the  q-  and  d-axis  currents  are  not  limited 
symmetrically.  Instead,  the  d-axis  current  is  given  priority  so 
as  to  maintain  the  desired  flux  level.  This  is  important 
because  it  prevents  a  transient  which  invokes  the  current 
limit  from  disturbing  the  rotor  flux  and  thereby  resulting  in 
loss  of  proper  flux  orientation. 

Based  upon  the  limited  current  command  and  the 
mechanical  rotor  speed  CO™,  the  indirect  field  oriented 
control  illustrated  in  Fig.  10  establishes  the  speed  and 
position  of  the  synchronous  reference  flame,  denoted  co«and 
0 , ,  respectively.  Note  that  part  of  the  indirect  field  oriented 
control  which  establishes  the  q-  and  d-axis  current 
commands  has  been  incorporated  into  the  link  stabilizing 
control. 


Figure  8.  Link  stabilizing  control. 

In  addition,  although,  the  details  will  not  be  set  forth 
herein,  the  fact  that  indirect  field  oriented  control  makes  use 
of  the  drive  parameters  coupled  with  the  size  of  the  drive 
neccsritatfM;  the  use  of  an  on  line-parameter  estimator  such 
as  a  Kalman  filter.  This  aspect  of  the  control  has  not  been 
set  forth  herein  but  is  an  important  part  of  the  indirect  field 
oriented  controller. 

The  final  aspect  of  the  control,  and  the  one  which  will  be 
the  focus  of  this  paper,  is  the  design  of  the  current  regulator. 
The  function  of  the  duty  current  regulator  will  be  to  establish 
the  duty  cycle  signal  to  a  sine-triangle  modulator  such  that 
the  actual  currents  are  equal  to  the  commanded  currents. 
Inputs  to  the  current  regulator  include  the  commanded  q- 
and  d-axis  currents,  the  dc  link  voltages,  the  speed  and 
position  of  the  synchronous  reference  flame,  and  finally  the 
instantaneous  phase  currents. 

VL  Decoupling  Transformation 

From  Fig.  4,  it  is  apparent  that  the  q-  and  d-axis  circuits 
are  highly  coupled  which  complicates  the  control  design. 
However,  by  suitable  selection  of  a  transformation,  the  q- 
and  d-axis  sets  can  be  decoupled.  In  order  to  derive  this 
transformation,  it  is  convenient  to  first  express  the  stator 


Figure  7.  Torque  limiter. 


Figure  9.  Current  limit  control. 


where 


Figure  10.  Field  oriented  controller. 

voltage  equations  in  subtransient  form.  In  the  case  of  the  q- 
axis,  solving  (24,26)  for  the  rotor  current  and  substitution 
into  (19)  yields 

Xqsi  —  L  Isi qsf  +  {L^lr  1 w)  Z  / qsf  + 

Similarly,  for  the  d-axis 

Xcisf  —L  isi dsi  +  {L^ir  I IZ  m)  Z  Z  dst  + 


By  manipulation  of  (12,13,17,18,28,29)  it  is  possible  to 
express  the  stator  voltage  equation  solely  in  terms  of  the 
stator  currents  and  rotor  flux  linkages.  In  particular, 

■r  •  T 

-rs\qs^^i  +<r>e^Z5»^iri23  ^  ^ 


U3.3i,. 


CD,(Ltlli„)U3.3i<i.,.  -%-U3.iV,,  +  co.^U3.iXrf.  (30) 
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and 
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U3.3^<^t23 


03.(i;ilI™)U3,3i,...  -^U3.iV,,-co.|^U3.iV  (31) 


and 
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1 
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(35) 

(36) 

(37) 


(38) 


Note  that  (38)  is  almost  orthogonal,  except  that  the  last  row 
has  been  modified  so  that  the  z-component  of  the  q-  and 
d-axis  variables  would  correspond  to  the  average  of  those 
variables.  This  transformation  yields 


''q^yyy  =  ^  Luj)i<p,y.  +  (OeLuids,yy+ 


J6 

0)«(Z»(rl 


fe)  ''r  +  (f-mli^Zr)^  €3.31^5^+ 

Lm)C3Mds^  “  qr-^^rJT^Xl'^dr 


(39) 


and 


yds^  =  “  (i^eLlsidS:^^ 
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^  6 


Lm)CxZ^qSjyt  ” 


(40) 


where  Cy  is  an  i  by  j  matrix  whose  i  j*th  element  is  one  and 
whose  remaining  elements  are  zero.  Upon  consideration  of 
the  definition  of  it  can  be  seen  the  x-,  y-,  and  z- 
components  of  the  q-  and  d-axis  stator  voltage  equations  are 
imcoupled. 


Vn.  Xy-Component  Control 


The  control  of  the  x-axis  component  of  the  q-  and  d-axis 
currents  is  as  foUows.  First,  extracting  the  x-axis 
components  from  (39)  and  (40), 


y qsx  ^ s^qst  '^^e^ls^dsx  (41) 

Vdsx  =  ^^^dsx  -(^eLlslqsx 


where  Ug  is  a  /  row  by  j  column  matrix  of  ones  and 


f<isi23  \_  fds\  fdS2  fdsi  J 


(32) 

(33) 


where  f  may  be  a  voltage,  current,  or  flux  linkage. 

In  order  to  decouple  the  three  sets  of  q-axis,  it  is 
convenient  to  transform  (30-3 1)  using  the  transformations 


^qsjtr  “ 


One  feature  of  interest  in  (41-42)  is  the  coupling  between  the 
q-  and  d-axis,  which  is  troublesome  because  it  introduces  a 
lightly  damped  eigenvalue  whose  natural  frequency  is  a 
function  of  operating  point  through  co  ^ .  A  suitable  control 
law  which  eliminate  this  problem  is  to  command  a  x- 
component  of  the  voltage  of 

=  (HeLisidsx  +  {Kpqx  +  qs^  ^  (^3) 

^dsx  ”  ^eLlsiqsx  {^pdx  +  “^)(J ds^  ~  ^dsj)  (^^) 


(34) 


where  s  is  the  Laplace  operator  and  Kpqx  Kjqx  denote  the 
x-component  of  the  q-axis  proportional  and  integral  gains, 
respectively.  This  control  igorithm  is  depicted  in  Fig.  11. 
Assuming  that  the  desired  voltage  is  obtained,  the  transfer 
function  between  the  q-axis  current  command  and  the  q-axis 
current  becomes 

_ KpqxS-^Kjqx _ 

i gjj,  '^K.pqs^S  +  K.icfx 

From  (45),  it  can  be  seen  that  the  poles  may  be  arbitrarily 
placed.  If  the  poles  are  widely  spaced,  and  the  fastest  pole  is 
faster  than  the  stator  leakage  dynamics,  the  zero  will  nearly 
cancel  the  first  pole,  so  that  the  second  pole  dominates  the 
response.  The  transfer  function  for  the  d-axis  is  identical, 
and  the  control  of  the  y-axis  component  is  identical  to  that  of 
the  x-axis  component  of  the  current.  Typically,  the  y-axis 
component  gains  will  be  identical  to  those  of  the  x-axis. 
However,  the  d-axis  gains  may  occasionally  be  selected  to  be 
different  than  that  of  the  q-axis. 


Vra.  Z-AXIS  CONTROL 


q- 


The  voltage  equations  governing  the  z-component  of  the 
and  d-axis  may  be  expressed  as 

T  ■  T 

^ ^^eLzlqSz 
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(46) 


di 


Vdsz  =  ^zidsz  -  <^eLzi(bz 

in  ^dr 
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"qr 


(47) 


^ds 


where 

r.  =  r,  +  76(fr)V,  (48) 

Lz  =  L,s+I^{L\,\\L„)  (49) 

The  dynamics  of  the  z-component  of  the  q-  and  d-axis 
currents  is  similar  to  the  design  of  the  dynamics  of  the  x- 
and  y-components  of  the  q-  and  d-axis  currents,  with  the 
exception  that  the  z-component  is  coupled  to  the  rotor 
circuits.  However,  the  stator  voltages  and  currents  vary  on  a 
time  scale  much  faster  than  that  of  the  rotor  flux  linkages. 
Therefore,  the  rotor  flux  linkage  terms  may  be  treated  as 
disturbances.  A  suitable  control  strategy  is  thus  to 
synthesize  the  commanded  z-components  of  the  voltages  as 

v;,,  +(^wr +  (50) 

+  {Kp^  +  ¥)(4.  -  '■*.)  (5 1) 

which  is  identical  to  the  form  of  the  control  for  the  x-  and 
y-components  of  the  currents.  Again  assuming  the  actual 
voltage  are  equal  to  the  commanded  voltages  the  transfer 
fimction  of  the  z-component  of  the  stator  current  may  be 
expressed 

hst  _  KpqzS'^Kjqz _  ^ 

^  qsx  LzS^  ^iqz 

As  with  the  case  of  the  x-  and  y-  component  controllers  the 
poles  may  be  arbitrarily  placed,  and  the  pole  placement  of 
the  q-  and  d-axis  may  be  placed  independently,  if  so  desired. 


DC  Zero  Sequence  Control 


Unlike  a  wye-connected  3 -phase  machine  zero  sequence 
currents  can  be  present  in  the  drive  system  shown  and  should 
be  actively  controlled  to  zero.  The  control  design  of  the  zero 
sequence  currents  is  straightforward  since  each  zero 
sequence  circuit  is  dynamically  equivalent  to  a  xmcoupled  RL 
circuit.  Herein,  a  simple  first  order  low  pass  controller  of  the 


form 


Vnji  — 


Ko 


XoS-^  1 


(53) 


is  utilized,  wherein  'n'  may  be  'O',  'a',  or  'P*.  Although  each 
zero  sequence  circuit  could  have  a  different  gain,  there  is  no 
reason  why  this  would  be  desirable  since  the  zero  sequence 
circuits  are  all  identical.  The  low-pass  filter  portion  of  the 
control  helps  filter  out  switching  noise.  It  has  been  found 
through  computer  simulation  that  the  zero-sequence  circuits 
are  much  more  susceptible  to  switching  noise  than  the  q-  and 
d-axis.  Based  on  (46),  the  characteristic  equation  of  the  zero 
sequence  current  is  given  by 


Figure  1 1 ,  QD  set  current  control. 


-  fLis^rsXo']  fr^Ko')  ^ 


(54) 


Figure  12.  Current  Regulator. 


Using  (53),  by  suitable  selection  of  to  and  Ko  the  poles  may 
be  arbitrarily  placed.  There  are  no  zeros  associated  with  the 
zero  sequence  dynamics, 

X.  Current  Regulator  Design 

The  finalized  current  regulator  is  illustrated  in  Fig.  12. 
Therein  the  limited  desired  qd-  currents  are  resolved  into 


x,y,  and  z-components.  The  measured  currents  are  also 
transformed  to  qd  variables,  and  then  divided  into  xyz 
components.  Based  on  the  speed  of  the  synchronous 
reference  frame,  the  commanded  xyz  currents  and  actual  xyz 
currents,  the  xyz  components  of  a  voltage  commanded  is 
determined,  from  which  the  qd  voltage  command  is  found. 
Finally,  the  qd  voltage  command,  along  with  the  zero 
sequence  voltage  commands  are  transformed  to  machine 
variables. 


T 


A  sine-triangle  modulator  is  used  to  achieve  the  desired 
voltage  command  for  each  phase.  In  particular,  a  duty  cycle 
is  formulated  for  each  phase  (denoted  di)  which  is  compared 
to  a  unity  amplitude  triangle  wave  (denoted  ij)  as  illustrated 
in  Fig.  13.  The  switching  signal  is  formulated  in  accordance 
with 


[  -1  if  di  <  ti 


(55) 


The  triangle  wave  signals  of  adjacent  phases  are  phase 
displaced  from  each  other  by  an  equal  amount  and  the 
frequency  of  the  triangle  waves  is  2  kHz.  It  can  be  readily 
shown  that  upon  neglecting  conduction  drops,  the  fast 
average  of  the  phase  voltage  is  given  by 

Vi  =  diVjej  (56) 

where  j  denotes  the  dc  link  corresponding  to  the  /'th  phase. 
Inverting  this  relationship  provides  a  method  of  achieving 
the  desired  voltage  as  depicted  in  Fig.  13. 

After  selection  of  suitable  pole  locations,  the  current 
regulator  gams  and  other  control  system  parameter  were 
determined  and  appear  in  Table  3 . 


Figure  13.  Sine-triangle  modulation  signals  and 


Table  3.  Induction  Motor  Controller  Parameters 

r^y.^  =  403  kNm 

r.../^=403MNin/s 

o 

o 

II 

j 

^./,>,  =  650A 

v^,  =  5618V 

r- 20  ms 

K=n.oivik 

pax 

K.  =  306.6  V/As 

tax 

=  61.65  V/A 

K  =  17.07  V/A 

pax 

=  306.6  V/As 

K.  =  504.9  V/As 
_iS _ 

A:^  =  4.175 

r^  =  4  ms 

X^*=53  Vs 

